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ABSTRACT 


The  results  of  a  five-month  analytical  program  are  presented^ 
describing  a  stress  analysis  method  for  a  three-dimensional 
fiber  composite  structure.  The  Intersections  of  the  elements 
or  webs  of  this  structure  form  polyhedra  in  the  form  of  tri¬ 
angular  prisms,  Tetra-core,  and  truncated  forms  of  these. 

The  finite-element  model  generates  the  elements  and  polyhedra 
automatically  with  an  optimization  routine  which  selects  the 
minimum  weight  structure  that  will  meet  the  stiffness  and  loads 
requirements. 
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INTRODUCTION 


This  docairant  contains  a  description  of  a  three-dimensional 
finite- oiom.ent  stress  analysis  computer  model  developed  to 
provide  a  tool  for  defining  the  geometry  of  Tetra-core 
structures.  The  model  uses  m.aterial  properties  and  applied 
loads  in  determining  the  final  composite  configuration.  The 
primary  polyhedron  shape  is  a  tetrahedron;  one  form  of 
the  model  can  produce  a  Tetra-core,  v.-here  all  elements  inter¬ 
sect  to  form  an  equilateral  tetrahedron. 

The  objective  of  the  program  was  to  develop  a  finite-element 
stress  model  that  would  provide  a  method  of  analysis  of  a 
compicx  nonhomoqenous  structure.  The  model  weis  developed 
to  run  on  tlie  IBM  360/44.  The  size  of  problem  that  can  be 
run  is  limited  only  by  the  number  and  size  of  high-speed 
discs  available. 

Tetra-core  is  a  new  construction  concept  developed  for  use  on 
helicopters,  airplanes,  and  other  vehicles  where  lightweight 
and  damage -tolerant  structures  are  desirable.  A  typical 
application  is  shown  in  Figure  1.  It  is  made  by  filament 
winding  fiber  glass  or  advanced  composite  fibers  into  an 
interlocking  group  of  tetrahedrons,  as  shovm  in  Figure  2. 

This  arrangement  appears  to  give  a  liqhtv/eight,  truss-type 
structure  that  will  be  able  to  sustain  damage  and  still  carry 
a  load,  because  of  the  high  degree  of  structural  indeterminancy. 
Tetra-core  is  constructed  using  a  standard  filament- winding 
machine.  The  fibers  are  laid  in  legs  oriented  in  one  direction, 
then  the  fibers  of  the  next  direction  are  laid  over  them,  thus 
providing  good  load  transfer  between  legs.  Adhesive  is  applied 
to  the  fiber  as  it  is  being  wound.  The  final  composite  is  then 
ready  for  the  curing  process. 

Several  stress  analysis  methods  have  been  considered  for  use 
v,’ith  Tetra-core.  A  simple  stress  analysis,  in  which  the  load 
in  each  leg  is  resolved  into  its  components  by  a  summation  of 
forces,  does  not  account  for  the  structural  indeterminancy  of 
Tetra-core,  and  will  probably  give  conservative  results. 

The  classical  laminated  plate  theory  widely  used  with  boron 
and  graphite -layered  panels  does  not  apply  to  Tetra-core  since 
it  is  not  a  homogenous  material.  Each  leg  of  a  Tetra-core  can 
be  analyzed  using  laminated  plate  theory,  but  a  method  that 
accounts  for  the  redundancy  caused  by  the  interv/eaving  of  the 
legs  should  be  used  in  determining  the  strength  of  a  Tetra-core 
panel . 
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The  finite-element  method  was  used  for  this  contract,  along 
with  an  automatic  grid  generator  to  automatically  idealize 
the  staructure  and  apply  leads  to  it. 

Two  types  of  elements  are  used  in  this  program:  a  linear 
strain  triangle  (6  nodes)  for  the  legs  of  the  Tetra-core,  and 
a  constant  strain  triangle  (3  nodes)  for  the  face  sheets 
attached  to  the  core.  The  linear  strain  triangle  is  required 
for  the  logs  because  of  the  change  of  strain  from  tension  to 
compression  between  top  and  bottom  of  the  tetrahedron  when 
the  panel  is  subjected  to  a  bending  moment.  If  a  constant 
strain  triangular  element  were  used  in  this  strain  field, the 
result  v;ould  be  an  averaging  of  the  maximum  and  minimum 
stresst’s  since  the  element  would  enforce  a  constant  strain 
over  the  triangle.  A  constant  strain  triangle  may  be  used 
for  the  face  sheets,  however,  since  the  stress  field  does  not 
contain  changes  of  sign  or  large  strain  gradients.  This  allows 
a  mere  efficient  use  of  the  available  core  storage  space  since 
additional  nodes  do  not  have  to  be  generated  at  the  midpoint 
of  each  side. 

The  program  has  been  v/ritten  to  automatically  generate  finite - 
element  models  for  Tetra-core  structures  made  up  of  flat 
plates,  cylinders,  or  airfoils.  T\v’o  basic  types  of  Tetra-core 
may  be  generated:  a  true  Tetra-core,  where  the  sides  of  the 
Tetra-core  meet  at  a  point  to  form  a  tetrahedron ;  and  a 
truncated  Tetra-core,  in  which  the  tops  have  been  "cut  off"  of 
the  tetrahedron.  The  basic  tetrahedron  geometry  can  be  varied 
to  give  a  skewed  tetrahedron.  Self-equilibrating  nodal  loads 
are  automatically  applied  to  the  model  based  on  the  input  com¬ 
bination  of  in-plane  and  out-of-plane  loads.  Multiple  load 
cases  may  be  run.  Boundary  conditions  (nodal  fixities)  are 
applied  to  the  model  to  hold  it  from  translating  or  rotating 
without  providing  reactions  for  the  nodal  loads.  This  com¬ 
bination  of  load.s  and  boundary  conditions  results  in  the 
idealization  of  a  coupon  of  Tetra-core  to  which  loads  corre¬ 
sponding  to  those  in  some  location  of  a  larger  structure  are 
applied.  A  conventional  finite-element  analysis  would  be  run 
of  the  total  structure  or  a  part  of  the  structure.  Stresses 
from  this  run  would  then  be  used  to  apply  loads  to  the  Tetra- 
core  model  corresponding  to  critical  locations  on  the  total 
structure. 

Deflections  are  calculated  using  a  Choleski  triangularization 
method  with  substitution  to  solve  the  equation  [K]  {5}  =  {P}. 

Stresses  and  margin  of  safety  are  calculated  at  several  loca¬ 
tions  in  the  linear  strain  elements.  A  buckling  analysis  of 
a  triangle  in  each  Jeg  is  done  using  the  buckling  equation  of 
an  orthotropic  rectangle  of  the  .same  size  as  a  triangle  in  that 
leg,  since  no  orthotropic  triangle  buckling  equation  are  avail¬ 
able.  An  optimization  of  the  Tetra-core  geometry  can  be  made 
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using  a  steepest  descent/side-step  method  to  find  a  minimum 
weight  design  for  the  given  loading  conditions.  A  nonlinear 
analysis  can  be  run  using  stress-strain  curves  of  the  material 
in  each  leg  to  account  for  the  change  in  stiffness  of  each 
plate  element  with  increasing  load.  The  effect  of  a  hole  in  a 
Tetra-core  structure  can  be  analyzed.  The  program  will  give  a 
zero  stiffness  to  plates  connecting  to  a  node  at  the  center  of 
the  hole.  This  causes  the  lo^^ds  to  be  redistributed  into  plates 
around  the  hole,  creating  the  effect  of  a  stress  concentration. 

The  program  is  now  configured  for  the  IBM  360/44,  but  it  can 
be  easily  modified  to  run  on  any  360.  Converting  to  a  CDC 
6600  would  be  more  difficult,  since  direct-access  disc  drives 
are  required. 

Figure  3  shows  the  elements  of  the  program  and  the  sequence  of 
calculation.  Each  program  element,  input  and  output  data,  and 
model  usage  are  described  in  the  following  chapters.  This 
includes  input  format,  output  format,  and  a  complete  listing  of 
the  program. 
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Figure  1  Sequence  of  Operetiora  in  Tetre-Core  Analysis  Progranv 
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MODEL  DESCRIPTION 


This  program  has  been  written  to  automatically  generate  a 
finite-element  model  as  a  series  of  node  points  in  space  -and 
connect  them  with  plate  elements  into  an  idealization  of  a 
Tetra-core  structure  (Figure  1) .  It  will  automatically  apply 
nodal  loads  to  the  model  to  represent  uniform  applied  loads. 
Tetra-core  flat  plates,  cylinders,  and  airfoils  can  be 
analyzed.  The  program  uses  the  finite-element  model  it  has 
generated  to  calculate  stresses  and  deflections  for  each 
element  and  each  node.  A  margin  of  safety  is  calculated  for 
each  plate  element.  A  simple  buckling  analysis  is  performed 
on  a  plate  in  each  leg.  An  optimization  can  be  done  to 
determine  the  least-weight  Tetra-core  design  for  a  given  set 
of  loads.  A  nonlinear  analysis  can  be  performed  by  inputting 
a  set  of  material  stress-strain  curves  and  applying  loads  in 
increments  to  compute  the  degradation  in  each  plate  element. 
The  effect  of  a  hole  in  the  model  can  be  simulated  by  giving 
a  low  stiffness  to  the  plates  connected  to  a  specified  node. 

A  diagram  of  the  sequence  of  computations  performed  by  the 
program  is  given  in  Figure  3. 

INPUTS 

Input  data  required  to  run  the  Tetra-core  analysis  program 
has  been  kept  to  a  minimum  to  simplify  the  vi;ork  required  for 
trade-off  studies  on  the  effect  of  geometry,  materials,  and 
other  variables.  Basic  inputs  are: 

o  Type  of  analysis  to  be  run — nonlinear,  optimization, 
damage,  face  sheets,  flat  plate,  cylinder,  airfoil,  and 
true  or  truncated  tetrahedron 

o  Optimization  controls 

o  Size  of  specimen  to  be  analyzed  and  tetrahedron  geometry 

o  Material  data — moduli  and  Poisson's  ratio,  allowable 
stresses,  and  stress-strain  curves 

o  Loading  to  be  imposed  on  the  model 

Using  these  data,  the  program  generates  a  set  of  node  points 
and  plates  to  model  the  Tetra-core  element  and  applies  the 
specified  loadings  as  nodal  loads. 

The  basic  tetrahedron  geometry,  shown  in  Figures  4  and  5,  is 
input  as  the  length  of  a  side,  the  actual  tetrahedron  height, 
the  theoretical  height  at  which  the  apexes  would  meet,  and 
the  offsets  which  determine  the  amount  of  deviation  from  a 
true  tetrahedron. 
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FRONT  VIEW 
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Figure  4.  Basic  Tetrahedron  Geometry. 
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If  the  actual  height  and  the  theoretical  height  are  identical, 
then  a  true  tetrahedron  will  result.  Two  generation  options 
are  available  in  the  program:  a  true  tetrahedron  and  a  trun¬ 
cated  tetrahedron.  The  true  tetrahedron  option  should  be  used 
whenever  possible,  since  it  requires  fewer  nodal  points  and 
plates  to  be  generated,  and  uses  a  more  efficient  nodal 
numbering  system  that  allows  a  smaller  bandwidth  when  used 
in  a  cylinder  or  an  airfoil.  If  the  theoretical  height  for  a 
truncated  tetrahedron  is  input  as  a  large  r:~ber  (such  as 
10  inches) ,  then  a  three-sided  figure  Vvith  near-vertical  walls 
will  result.  The  truncated  Tetra-core  model  cannot  be  used 
to  generate  a  true  Tetra-core,  since  three  node  points  at  the 
top  of  a  tetrahedron  would  have  the  same  cooidinates .  This 
would  cause  a  singular  stiffness  matrix  and  result  in  a  bad 
solution.  Offsets  can  be  used  to  change  the  shape  of  the 
figure  from  a  true  tetrahedron  v/ith  angles  of  60  degrees,  if 
no  offsets  are  used,  to  an  isosceles  tetrahedron  of  any  angle, 
if  offset  is  used,  to  a  skewed  tetrahedron,  if  0  is  used. 

The  number  of  legs  on  the  X  and  Y  faces  (LX  and  LY) must  be 
even  numbers  (2,  4,  6,  etc.). 

Plates  to  connect  the  node  points  and  form  a  Tetra-core  model 
are  generated  sc  that  the  local  X-axis  is  always  in  the 
direction  of  the  fibers  in  that  leg,  as  shown  in  Figure  6. 

For  sandwich  faces,  the  local  X-axis  is  in  the  global  X  direc¬ 
tion.  Elastic  moduli,  allowable  stresses,  and  stress-strain 
curves  are  input  separately  for  each  leg  and  each  face  sheet, 
thus  allowing  each  leg  to  be  a  different  material  if  desired. 
Material  stress-strain  curves  are  input  only  if  a  nonlinear 
analysis  is  run. 

Using  the  input  data  LX  and  LY  and  the  geometry,  the  tetra¬ 
hedrons  are  automatically  integrated  into  a  Tetra-core  node 
system  as  shown  in  Figure  7  for  a  true  tetrahedron  and  in 
Figure  8  for  a  truncated  tetrahedron.  The  node  point  numbering 
system  is  shown  for  each  type.  Other  examples  are  shown  in 
Appendix  I. 

The  number  of  plates  generated  depends  on  the  input  number  of 
legs.  Plates  are  divided  into  Vertical,  Skev;  A, and  Skew  B 
legs  by  the  program.  In  Figures  9  and  10,  the  Skew  A  and  B 
plates  generated  for  a  true  tetrahedron  with  10  legs  on  the  X 
face  and  6  legs  on  the  Y  face  are  shown.  Plates  are  numbered 
in  sequence  with  plates  in  the  Vertical  legs  generated  first, 
then  plates  in  Skew  A  legs,  then  plates  in  Skew  B  legs.  t-Jhen 
cylinders  and  airfoils  are  generated,  several  additional  plates 
are  added  to  "zip  up"  the  seam  and  connect  the  two  sides  v/ith 
more  Skew  A  and  Skew  B  plates.  Figure  11  shows  the  Skew  A 
plates  making  up  a  cylinder.  Face  sheets  can  be  automatically 
added  to  the  true  Tetra-core  model  to  give  a  sandwich  element. 
It  was  not  considered  feasible  to  generate  a  face  sheet  for 
the  truncated  Tetra-core  element  due  to  the  complex  shapes 
which  would  be  required. 
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VERTICAL  LEG 


Local  X-Y  Axis  of  Platn  li  In  the 
Plane  of  the  Leg  With  the  X  Axle 
in  the  Fiber  Direction  and  the  Y 
Axis  at  90P  in  tha  Plane  of  the 
Plate 


Moduli  and  Allowible  Strews 
Are  Input  for  Each  Lag  in  the 
Local  X'Y  Coordinate  System 


Figure  6.  Material  Properties  Required  at  Input. 


•  Nodes  Connected  by  Dark  Lines  Are  on 
Upper  Surface  of  Plate 

•  Nodes  Connected  by  Light  Lines  Are  on 
Lower  Surface  of  Plate 


1 

X 

Figur«7.  Nodal  NumbaringSyitwn,  True Tetra-Core  Fist  Plata,  =  4,  Ly  =  4. 
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FigurtS.  Nodal  Numbering  Syttwn.Truncatid  Tvtrs-Coft  Flat  PIfti,  *  4,  Ly  4. 
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firm  Skaw  A  PM* 
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First  Plate  Element 
io  Skew  B  Lagi 


Figure  11.  Skew  A  Plates  in  Cylinder  (Computer  Plot) . 
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When  generating  a  cylinder,  the  program  transforms  the  flat 
plate  Y  and  Z  coordinates  into  a  circular  shape  with  a  circum¬ 
ference  equal  to  the  Y  dimension  of  the  plate.  Thus  the  out¬ 
side  radius  of  the  cylinder  is  the  X  dimension  of  the  plate 
divided  by  Ztt.  The  airfoil  section  is  generated  with  an 
input  set  of  thickness/chord  and  X/chord  ratios  to  define  its 
shape.  Vertical  legs  are  defined  to  be  parallel  to  the  X-axis 
for  cylinders  and  airfoils  to  give  the  maximum  bending  stiff¬ 
ness  . 

Node  point  loads  are  automatically  applied  to  the  finite- 
element  model  based  on  the  input  loading  conditions.  Loads 
are  automatically  applied  in  self-equilibrating  sets;  i.e., 
a  load  applied  to  one  side  is  reacted  by  an  equal  and  opposite 
load  on  the  opposite  side.  Loading  of  the  linear  strain  tri¬ 
angle,  which  is  used  for  the  Vertical,  Skew  A,  and  Skew  B  legs, 

includes  a  load  on  the  mid-node  point  between  the  vertex  node 
points.  This  load  is  four  times  the  load  at  a  vertex  for  a 
uniform  applied  load.  The  vertex  node  point  loadings  generated 
by  inputting  flat  plate  end  loads  in  the  X  and  Y  directions  and 
shear  loads  for  the  true  tetrahedron  case  are  shown  in  Figures 

12a,  b,  and  c.  fbr  moment  loads,  the  load  on  the  center  node 

point  of  the  triangle  side  is  zero,  since  the  two  vertex  loads 
cancel  out.  Node  point  loadings  generated  from  input  moments 
in  the  X,  Y,  and  XY  directions  are  shown  in  Figures  12d,  e, 
and  f  for  the  true  tetrahedron.  Out-of-plane  shear  loads  can 
be  applied  to  the  model.  A  moment  is  automatically  applied 
to  react  the  moment  caused  by  the  shear  loads  as  shown  in 
Figures  12g  and  h  for  vertical  shear  loads  on  the  X  and  Y  faces. 

The  truncated  tetrahedron  is  loaded  in  a  similar  manner,  except 
that  the  lower  node  points  in  an  edge  are  given  twice  the  load 
of  the  upper  node  points,  since  there  are  twice  as  many  upper 
node  points.  The  mid-side  node  point  load  is  four  times  the 
upper  node  point  load.  More  detailed  examoies  of  loads  applied 
to  Tetra-core  flat  plates  are  shown  in  Appendix  I. 

The  cylinder  model  can  be  loaded  with  moment,  torque,  and  end 
loads.  The  input  load  is  used  as  the  total  load  on  the  section 
and  is  divided  among  all  nodes  on  that  end.  The  nodal  loads 
generated  from  an  input  torsion  load  on  a  cylinder  are  shown 
in  Figure  13. 

The  airfoil  model  can  be  loaded  with  moment,  torque,  and  end 
loads  in  a  similar  manner  to  the  cylinder.  A  torsion  load 
applied  to  a  diamond  airfoil  section  is  shown  in  Figure  14. 
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Figure  1Z  Flat  Plate  Applied  Loads  . 
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(e)  Y  MOMENT,  PMY 


(gl  X  SHEAR,  XQSHR 


(h)  Y  SHEAR,  YQSHR 


Figure  12.  Continued. 
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Figure  14.  True  Ti 


Loads  can  be  applied  to  the  true  Tetra-core  specimen  being 
analyzed  in  any  combination  (end  loads,  moments,  vertical 
shears,  etc.).  The  truncated  Tetra-core  specimen  can  be 
loaded  with  end  loads  and  moments  in  any  combination,  but  if 
a  vertical  shear  load  is  applied,  then  only  certain  combina¬ 
tions  of  loads  may  be  applied  at  the  same  time,  because 
of  the  fixities  that  are  applied  to  the  four  sides  of 
the  panel  to  keep  the  model  stable.  If  no  vertical  shear 
loads  are  applied,  the  Z-freedcn  is  fixed  for  all  nodes  on 
the  panel  sides.  If  a  vertical  shear  is  applied  on  the  X 
face,  the  side  nodes  are  fixed  in  the  Y  direction,  which 
permits  movement  in  the  Z  direction.  A  vertical  shear  load 
on  the  Y  face  causes  edge  nodes  to  be  fixed  in  the  X  direction. 
Thus,  if  a  load  in  the  Y  direction  is  applied  to  a  panel  with 
a  vertical  shear  load  on  the  X  face,  the  effect  of  the  Y  loads 
will  be  lost,  since  they  are  applied  to  nodes  that  cannot  move 
in  the  direction  of  the  load. 

Moment,  torsion,  and  end  loads  can  be  applied  in  any  combination 
to  both  true  and  truncated  Tetra-core  cylinders  and  airfoils. 

The  effect  of  a  hole  in  the  Tetra-core  can  be  simulated  by 
giving  the  plates  connected  to  a  node  in  the  center  of  the 
panel  a  zero  stiffness.  This  unloads  those  plates  and  transfers 
the  load  to  the  remaining  plates,  providing  the  effect  of  a 
stress  concentration  around  the  hole. 

AUTOMATIC  MODEL  GENERATOR 

The  Tetra-core  analysis  program  has  been  written  so  that  a 
finite-element  model  will  be  generated  from  the  input  data. 

o  Type  of  Analysis — Card  Type  1 

o  Optimization  Controls — Card  Type  2 

o  Geometry — Card  Types  5,  6,  and  7 

o  Material  Data — Card  Types  4,  8,  9,  and  10 

o  Loading — Card  Types  12  and  13 

Flat  Plate  Nodal  Generator 

A  different  nodal  generator  subroutine  has  been  written  for 
each  type  of  Tetra-core,  true  and  truncated.  Subroutines 
NODGEN  and  NODGAN  are  used  to  generate  nodal  coordinates  for 
the  true  and  truncated  Tetra-core  models,  respectively.  True 
Tetra-core  nodes  are  generated  starting  at  the  Y-axis,  with 
numbers  increasing  as  the  X  distance  increases,  as  shown  in 
Figure  7.  Nodal  coordinates  (X,  Y,  Z)  on  the  bottom  surface 
of  the  Tetra-core  at  a  given  X  distance  are  generated  first. 
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using  the  repeating  geometry  of  Tetra-core,  then  those  on  the 
upper  surface.  Then  the  X  distance  is  increased  and  another 
cycle  of  numbering  is  begun.  For  the  truncated  Tetra-core,  a 
different  system  is  used  because  of  its  more  difficult  geom¬ 
etry.  Numbering  is  started  at  the  Y-axis  on  the  first  vertical 
leg,  as  shown  in  Figure  8.  Node  numbering  follows  the  first 
vertical  leg  to  its  end,  then  switches  to  the  second  vertical 
leg,  etc. 

Nodal  coordinates  are  generated  in  order  of  their  appearance 
on  the  leg,  resulting  in  a  mixing  of  nodes  from  the  upper  and 
lower  surfaces.  Nodes  in  the  last  vertical  leg  are  generated 
in  the  same  manner,  but  every  fourth  node  is  unconnected  to 
any  other  node  in  the  flat  plate  and  must  be  given  fixity  in 
the  X,  Y,  Z  directions.  If  a  cylinder  or  airfoil  is  generated, 
these  extra  nodes  are  used  to  connect  the  two  sides  of  the 
plate  after  it  has  been  folded  so  that  the  two  edges  almost 
meet. 


Cylinder  Nodal  Generator 


Cylinders  are  generated  from  the  flat  plate  coordinates 
previously  generated.  The  cylinder  is  assumed  to  have  a 
circumference  equal  to  the  length  of  the  X  face  of  the  flat 
plate  plus  the  width  of  one  tetrahedron.  This  gives  the 
cylinder  an  outside  radius  that  is  the  circumference  divided 
by  2it.  The  new  Y  and  Z  coordinates  of  each  node  are  calculated: 


6^  =  y^  flat  plate 


*27r/  circumference 


y^yj^  =  radius  *  cos  (6^) 

z  ,  =  radius  *  sin  (6 . ) 
cyl  3 


The  flat  plate  is  essentially  rolled  into  a  cylinder  so  that 
the  two  edges  of  the  flat  plate  are  separated  by  the  width 
of  the  one  tetrahedron  side  in  the  Y  direction.  The  same 
method  is  used  for  both  true  Tetra-core  and  truncated  Tetra- 
core.  An  example  of  the  end  of  a  true  Tetra-core  cylinder  is 
shown  in  Figure  13. 


Virfoil  Nodal  Generator 


Airfoil  nodal  coordinates  are  generated  by  transformation  of 
flat  plate  coordinates  using  a  curve-fitting  subroutine.  The 
same  basic  method  is  used  for  both  true  Tetra-core  and  trun¬ 
cated  Tetra-core.  The  flat  plate  already  generated  is 
essentially  folded  over  to  make  the  upper  and  lower  surfaces 
of  the  airfoil.  The  chord  length  of  the  airfoil  is  half  the 
length  of  the  X  face  of  the  original  flat  plate.  The  Z  co¬ 
ordinate  of  each  node  is  determined  using  the  curve- fitting 
subroutine  with  the  input  t/C  and  X/C  values.  The  distance 
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of  a  node  from  the  X-iixis  and  the  airfoil  chord  length  is  used 
to  calculate  the  X/chord  ratio  of  that  node.  This  X/chord 
ratio  is  used  in  the  curve-fitting  subroutine  to  determine  the 
thickness/ciicrd  ratio  of  a  point  on  the  surface  of  the  airfoil 
from  the  iiiput  t/C  ratios  of  the  airfoil.  The  Z  coordinate 
of  the  nodtj  is  tlion  the  tbickness/chord  ratio  for  that  point 
times  the  chord  length.  See  Figure  14  for  an  example  of  the 
end  of  a  diamond -shaped  airfoil  generated  by  the  program. 

Plate  Element  Generator 

Plate  elements  are  generated  automatically  to  connect  the 
nodes  from  the  nodal  generator  and  form  either  a  true  Tetra- 
core  or  a  truncated  Tetra-core  plate.  True  Tetra-core  plate 
elements  are  generated  in  subroutine  PLATGd.-  with  truncated 
Tetra-core  plate  elements  generated  in  subroutine  PLATGA.  The 
same  method  is  used  for  both  models,  but  a  different  subroutine 
is  needed  for  each,  since  tne  nodal  numbering  sequence  is 
different  for  each. 

Flal  Plate 

Plate  element  generation  starts  with  the  upper  plate  in  the 
left  vertical  leg.  Linear  strain  triangular  and  quadrilateral 
plates  are  used  for  the  Vertical  and  Skew  A  and  B  legs,  since 
they  can  be  expected  to  encounter  compression  strain  on  one 
side  of  the  leg  and  tension  strain  on  the  other  for  some  load 
cases.  The  mid-side  node  on  the  outer  face  of  the  element  is 
reduced  out,  since  it  will  not  be  connected  to  another  element 
and  would  cause  the  merged  stiffness  matrix  to  be  singular  if 
left  in.  Tv;o  linear  strain  triangles  are  added  to  get  the 
quadrilateral  linear  strain  element  used  in  the  truncated  Tetra 
core  model.  See  Figure  24  for  the  tria'igle  and  quadrilateral 
element  node  numbering  sequence. 

Plates  are  numbered  so  that  the  local  X-axis  of  each  plate 
element  is  in  the  fiber  direction  of  that  leg.  ’'Vertical  plates 
are  generated  by  following  the  vertical  leg  from  top  to  bottom, 
then  moving  to  the  next  vertical  leg  to  the  right,  and  starting 
at  the  top  again.  Each  vortical  leg  has  the  same  number  of 
vertical  elements. 

Skew  plate  elements  are.  then  generated,  starting  with  the  plate 
in  the  bop  left  of  the  flat  plate  segment,  as  shown  in  Figure  9 
Skew  A  legs  are  generated  from  upper  left  to  lower  right. 
Generally,  each  Skew  A  leg  has  a  different  number  of  plate 
elements.  Next,  Skev;  B  linear  strain  plate  elements  are  gene¬ 
rated,  starting  at  the  upper  right  corner  of  the  segment,  as 
shown  in  Figure  10.  Leg  generation  proceeds  from  upper  right 
to  lower  left,  with  a  different  number  of  plate  elements  in 
each  leg. 
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If  the  face  sheet  option  has  been  specified  ("IFACE"  =  1  on 
Card  1) ,  an  upper  and  lower  face  sheet  will  be  generated  for 
the  segment. 

Cylinder  and  Airfoil 

A  cylinder  is  generated  by  folding  a  flat  plate  model  into  a 
circular  cylindrical  segment.  The  two  edges  of  the  flat  plate 
facing  in  the  Y  direction  end  up  facing  each  other,  separated 
by  the  width  of  a  tetrahedron  side.  The  result  is  a  cylinder 
with  a  strip  missing  down  the  side.  Plates  are  added  to  connect 
the  two  sides  of  the  cylinder,  resulting  in  a  complete  cylinder, 
shown  in  Figure  11.  The  plates  making  up  a  Tetra-core  airfoil 
are  generated  in  the  same  way,  with  the  flat  plate  model  being 
folded  over  into  the  desired  shape  and  plate  elements  added  to 
"zip  up"  the  seam  where  the  two  ends  meet. 

Flat  Plate  Loads 


The  load  generator  takes  input  loads,  which  are  uniform  loads 
(given  in  either  pounds/inch  or  inch-pounds/inch) ,  and  applies 
them  as  nodal  loads  to  give  an  equivalent  total  load.  The 
same  method  is  used  for  both  true  Tetra-core  panels  and  trun¬ 
cated  Tetra-core  panels.  Since  a  different  nodal  numbering 
system  was  used  for  each  type  of  Tetra-core,  a  different  sub¬ 
routine  was  used  for  each.  Subroutine  LODGEN  generates  loads 
for  the  true  Tetra-core  model,  while  loads  for  the  truncated 
Tetra-core  model  are  generated  in  subroutine  LODGAN.  When 
generating  nodal  loads  for  the  linear  strain  triangle,  the 
mid-side  node  must  be  loaded  to  give  an  equivalent  uniform  load 
across  the  element.  Theoretically,  the  load  should  be  divided 
so  that  one-sixth  of  the  total  load  on  that  plate  goes  into 
each  corner  node  and  two-thirds  goes  into  the  mid- side  node^. 
This  is  done  automatically  in  the  program. 


In  the  case  of  a  uniform  load  in  the  X  direction,  the  input 
load  (pounds/inch)  is  multiplied  by  the  panel  width  and  divided 
by  the  number  of  nodes  on  that  face,  taking  account  of  the  fact 
that  some  nodes  have  more  load  than  others. 


PX  =  XLOAD  X  XLGT 
NXND 

where  XLOAD  =  Input  load  on  the  X  face  (Ib/in.) 

XLGT  =  Length  of  X  face 

NXND  =  Number  of  nodes  on  X  face 

NXND  =  2  X  LX  +  (2  X  LX  -  1)  x  4  for  true  Tetra- 
core  model 
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This  gives  the  X  direction  load  to  be  applied  to  each  node 
(PX) .  Loads  are  applied  to  the  apex  nodes  only  in  the  LODGEN 
and  LODGM  subroutines  since  mid-side  nodes  are  generated  in 
the  stress  analysis  overlay.  Nodal  loads  applied  to  give  the 
effect  of  a  uniform  load  in  the  X  direction  are  shown  in 
Figure  12a  for  a  true  Tetra-core  flat  plate. 

Figure  12b  shows  a  load  in  the  Y  direction  applied  to  the 
true  Tetra-core  panel.  The  nodal  load  PY  is  calculated: 


_  YLOAD  X  YL3T 
NYND 

where  YLOAD  =  input  load  on  the  Y  face  (Ib/in.) 

YLGT  =  Length  Y  face  (in.) 

NYND  =  2  X  LY  +  (2  x  LY  -1)  X  4  for  a  true  Tetra- 
core  model 

LY  =  Number  of  legs  on  Y  face 

For  a  uniform  in-plane  shear  load,  the  loads  shown  in  Figure 
12c  are  applied.  The  nodal  loads  are  calculated  in  the  same 
way  as  for  a  uniform  load  in  the  Y  direction,  except  that  a 
different  nodal  load  (PXYA  or  PXYB)  is  used  for  the  X  and  Y 
faces. 


PXYA  = 


XYLOD  X  XLGT 
NXND 


where  XYLOD  =  input  in-plane  shear  load  (Ib/in.) 

A  moment  load  in  the  X  direction  is  applied, as  shown  in 
Figure  12d,  for  a  true  Tetra-core  flat  plate.  For  a  moment 
load,  the  mid-side  nodes  are  not  loaded.  Nodal  loads  are 
calculated  using  the  moment  load  per  inch,  panel  width,  panel 
depth,  and  number  of  nodes  on  the  side  being  loaded. 


XMOM  X  J:LGT 
'  -HT-'x  LX - 

where  XMOM  =  Input  moment  on  X  face  (in. -Ib/in.) 

A  moment  load  is  applied  in  the  Y  direction,  as  shown  in 
Figure  12e.  The  nodal  loads  are  calculated: 


PMY 


YMOM  X  YLGT 
IIT  X  LY 


where  YMOM  =  Input  moment  on  Y  face  (in. -Ib/in.) 
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A  twisting  shear  moment  is  applied  as  shown  in  Figure  12f. 
Nodal  loads  PMXYA  and  PMXYB  are  calculated  separately  for  the 
X  and  Y  faces  of  the  panel: 


PMXYA  = 


XYMOM  X  YLGT 
HT  X  LY 


_  XYMOM  X  XLGT 
PMXYB  -  hT  X  LX 

where  XYMOM  =  Input  twisting  moment  (in.-lb/in.) 

Out-of-plane  shear  loads  can  be  applied  to  flat  panels  using 
the  input  loads  XQSHR  and  YQSHR.  A  combination  of  a  vertical 
shear  load  and  a  moment  load  is  applied  to  each  loaded  end  of 
the  panel,  as  shown  in  Figures  12g  and  12h.  The  moment  load 
PMX  is  applied  to  balance  the  moment  caused  by  the  two  vertical 
shear  loads  PZ.  The  resulting  set  of  loads  is  self  equili¬ 
brating,  so  no  reactions  will  result  at  the  fixed  nodes. 

Calculations  of  nodal  loads  for  a  given  X  and  Y  vertical  shear 
are: 


XQSHR  X  XLGT 

PZ  =  - - - 


_  XQSHR  X  SIDE  x  (LY  -  1)  X  XLGT 
^  "  2  X  HT  X  LX 

where  XQSHR  =  Input  vertical  shear  load  on  X  face  (Ib/in.) 


YQSHR  X  YLGT 
NYND 


_  YQSHR  X  SIDE  x  (0.86603  +  YZ  OFF)  x  (LX-1)  x  YLGT 

2  X  HT  X  LY 

where  YQSHR  =  Input  vertical  shear  load  on  Y  face  (Ib/in.) 

For  a  true  Tetra-core  flat  plate,  these  vertical  shear  loads 
can  be  applied  with  any  other  loads,  such  as  moment  loads  and 
in-plane  shear  loads;  etc.  for  the  truncated  Tetra-core  flat 
panel,  all  loads  cannot  be  applied  simultaneously.  A  vertical 
shear  load  YQSHR  can  be  applied  with  loads  and  moments  on  the 
X  face»  but  cannot  be  used  with  in-plane  shears  or  tv;isting 
moments  or  with  a  moment  on  the  Y  face.  Similarly,  the  verti¬ 
cal  shear  load  YQSHR  can  be  applied  v;ith  loads  and  moments  on 
the  Y  face,  but  not  with  shear  loads  or  moment  loads  on  the 
X  face.  The  reason  is  that  the  truncated  Tetra-core  model  is 


27 


unstable  on  the  edqes  and  must  be  supported  to  prevent  rigid 
body  rotation  in  some  of  the  plates.  This  is  done  by  fixing 
nodes  on  the  flat  plate  edges  in  the  Z  direction.  This  allows 
the  flat  plate  to  deflect  in  the  X-Y  plane  but  not  in  the  Z 
direction.  With  this  fixity,  no  vertical  shear  load  can  be 
applied  to  the  model  because  nodes  that  would  be  loaded  in  the 
Z  direction  are  fixed  in  that  direction  and  cannot  be  loaded. 
When  a  vertical  shear  load  on  the  X  face  is  specified,  the 
program  will  fix  the  cage  nodes  in  the  Y  direction,  instead 
of  the  Z  direction,  v/hich  allo’.vs  Z  nodal  loads  to  be  applied. 
Hov'ever,  loads  can  no  longer  be  applied  to  these  nodes  in  the 
Y  direction.  Therefore,  sivear  loads  and  mcmenc.  on  the  Y  face 
cannot  be  applied  at  the  same  time  that  vertical  shear  loads 
are  applied  to  the  X  face.  Similarly,  if  a  vertical  shear 
load  is  applied  to  the  Y  face,  all  edge  nodes  are  fixed  in 
the  X  direction,  and  shears  and  moments  on  the  X  face  cannot 
be  applied.  For  the  same  reason,  a  twisting  shear  moment  can¬ 
not  be  applied  in  conjunction  with  an  out-of-plane  shear  load. 

Cylinder  Loads  Generation 


Total  end  loads.  Lending  moments,  and  torsion  ioaos  are  applied 
to  cylinders  as  noflal  loads.  End  loads  are  applied  to  the 
cylinder  by  dividing  the  input  X  load  (in  pounds)  by  the  number 
of  nodes  on  one  end  of  the  cylinder,  taking  into  account  the 
fact  that  some  nodes  ha/e  a  larger  load  than  others. 


px  = 

MXN’D 

where  XLOAD  =  input  X  total  and  load  on  cylinder  (lb) 

A  torsion  load  is  applied  to  the  cylinder  as  nodal  loads 
tangent  to  the  surface  of  the  cylinder  at  the  node  being 
loaded,  as  shown  in  Figure  13.  Loads  on  the  inner  .surface 
of  the  cylinder  have  a  smaller  load  than  t)iose  on  the  outer 
surface  Ly  the  ratio  of  their  radii. 


PA 


TORQ _ 

LX'x  (r  ^ 
A 


X 


PB 


r 

"A 


where  TORQ  =  Input  total  torsion  load  on  cylinder  (in. -lb) 
r.  =  Radius  of  outer  surface  of  cylinder 
=  Radius  of  inner  surface  of  cylinder 
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A  bending  moment  is  applied  to  the  cylinder  as  a  tension  load 
on  nodes  below  the  Y-axis  of  the  cylinder  and  a  compression 
load  on  nodes  above  the  Y-axis  of  the  cylinder.  The  magnitude 
of  the  load  on  each  node  depends  on  the  distance  of  that  node 
from  the  Y-axis.  That  is,  the  neutral  axis  of  the  cylinder 
is  assumed  to  be  the  Y-axis,  and  a  straight-line  variation  of 
load  with  distance  from  the  neutral  axis  is  used. 


PMX 


XMOM 

sum' 


where  XMOM  =  Input  total  moment  on  end  of  cylinder  (in. -lb) 

nodes- 

SUM  =  4i:2  (Y.  ) 

j  =  l  ^ 


Airfoil  Loads  Generation 


End  loads  and  bending  moments  are  applied  to  the  airfoil  in 
the  same  way  they  were  applied  to  the  cylinder.  A  torsion 
load  is  applied  to  the  airfoil  by  applying  a  load  to  each  end 
node,  in  the  direction  of  the  airfoil  surface  at  that  node, 
as  shown  in  Figure  14.  Since  no  prediction  could  be  made  of 
the  location  of  the  shear  center  of  the  airfoil,  it  was 
necessary  to  apply  the  same  total  load  to  each  node.  This  was 
done  instead  of  ratioing  the  magnitude  of  each  load  by  its 
distance  from  the  shear  center,  as  was  done  for  the  cylinder. 

Nodal  Fixity 

In  a  finite -element  analysis,  it  is  necessary  to  restrain  the 
model  from  translating  or  rotating  as  a  rigid  body,  as  this 
will  result  in  a  singular  stiffness  matrix.  This  is  done  by 
removing  the  ability  of  some  nodes  to  move  in  the  X,  Y,  or  Z 
direction,  thus  "fixing"  them  in  that  direction. 

For  the  Tetra-core  flat  plate  analysis, nodes  on  the  Y  face 
next  to  the  X-axis  (nodes  3,  5,  11,  13,...,  29  in  Figure  7) 
are  fixed  in  the  Y  direction.  The  first  two  nodes  on  the 
first  vertical  leg  (nodes  3  and  5  in  Figure  7)  are  fixed  in 
the  X  direction.  The  lower  of  the  two  nodes  is  fixed  in  the 
Z  direction.  This  keeps  the  model  from  translating  or 
rotating  as  a  rigid  body. 

The  boundary  conditions  and  system  of  applied  loads  used  for 
the  flat  plate  Tetra-core  model  result  in  a  coupon  type 
analysis.  That  is,  the  plate  segment  being  analyzed  can  be 
thought  of  as  a  small  coupon  cut  out  of  the  larger  structure, 
with  loads  corresponding  to  the  loads  in  the  larger  structure 
at  that  point  applied  to  the  coupon  sides.  This  allows  the 
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use  of  a  coarse-grid  finite-element  analysis,  using  a 
standard  finite-element  program,  to  determine  the  loads  at 
any  point  in  the  structure.  Tetra-core  analyses  are  then 
run  of  several  points  in  the  structure ,  using  loads  found 
from  the  coarse-grid  analysis. 

For  the  truncated  Tetra-core  model,  additional  fixities  are 
applied  to  make  the  model  stable,  as  discussed  in  the  previous 
section.  Nodes  on  the  panel  edqe  are  given  fixity  in  the  Z 
direction  to  restrain  rigid-body  rotation  of  some  edge  plates. 
This  can  be  illustrated  by  looking  at  the  Skew  D  leg  in  the 
upper  right  corner  of  a  panel  (Figure  8)  .  It  can  be  seen  that 
this  plate  element  can  pivot  about  its  lower  surface  node 
(No.  46)  without  inducing  any  in-plane  strains  into  the  Skew  A 
or  vertical  plates  supporting  it.  This  results  in  a  singular 
stiffness  matrix,  but  can  be  prevented  by  restraining  movement 
of  the  edge  nodes  in  the  X,  Y,  or  Z  direction.  If  a  vertical 
shear  load  is  applied  on  the  X  faces  of  the  panel,  then  edge 
nodes  are  given  a  fixity  in  the  Y  direction  instead  of  the  Z 
direction.  If  a  shear  load  is  applied  on  the  Y  faces,  then 
edge  nodes  are  given  a  fixity  in  the  X  direction,  as  discussed 
in  the  Load  section. 

For  Tetra-core  cylinders  and  airfoils,  nodes  on  the  end  next 
to  the  X-axis  are  given  fixity  in  the  X,  Y,  and  Z  directions, 
providing  the  effect  of  a  cylinder  with  one  end  fixed  into  a 
stiff  foundation. 
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FINITE-ELEMENT  STRESS  ANALYSIS 


Linear  Strain  Element 


To  best  derive  the  stiffness  matrix  for  a  linear  strain,  six- 
node  triangular  element,  we  first  require  the  geometry  expressed 
in  terms  of  some  local  in-plane  coordinate  system. 


If  Xqp  =  X(Q)  - 

Yqp  =  Y(Q)  -  Y(P) 

Zqp  =  Z(Q)  -  Z(P) 


where  is  the  length  of  tHe>  line  PQ. 
Then  the  direction  cosines^  of  PQ  are 

h  ^  ^Qp/^1  >  0 

^2  ""  ^Qp/^'l 


^3  “  ^Qp/^1 

or  =  ^2  ~  ^3  ^1  5 


The  direction  cosines  of  the  normal  to  PQ,  through  R  are 
obtained  as  follows; 
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with  RR  =  Xj^p  +  Ypp  +  Xj 


and  the  direction  cosines  of  the  normal  are 


^£OC 


Referring  the  triangle  to  local  coordinates  shown  previously 
and  denoting  X,  „  -  X,  by  X.  ,  we  have 


a 


3 


^21  “  ^QP^l  ^QP^2 
^21  “  ^QP^l  ^QP^‘2 
^32  "  ^RQ''l  ■*■  ^RQ^2 
^32  ^RQ'^1  ^  '^RQ^2 
^31  ^RP^l  ^RP^2 


^31  ^RP'‘l  ■*■  ^RP^2 


+  Z 

+  z 

+  z 

+  z 

+  z 

+  z 


QP^3 

QP*‘3 

RQ^3 

RQ^3 

RP^3 

RP''^3 


2 

Following  the  method  of  Felippa,  the  stiffness  matrix 
of  the  six-node  linear  strain  triangle  may  be  expressed  as 


iiQ 

Ci2Q 

^13' 

2iQ 

C22Q 

^23' 

3iQ 

C32Q 

C33' 

where  B  relates  the  strain  at  any  point  within  the  triangle 
to  the  nodal  deflections, 
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and  Q  = 


2 


1 


1 


(3  X  3) 


12  1 
112 


^11 

S2 

^13 

f 

^1 

"2 

> 

^21 

^22 

^23 

i 

^2 

"12 

^  d 

Si 

S2 

^33 

“12 

^  d 

Using  Hooke's  Law  to  relate  stress  and  strain,  C23  = 
C31  =  Ci3  =  0,  and  further,  Cij  is  symmetric. 


Matrix  B  may  be 

expressed 

as  follows: 

3bi  -h^ 

4b2 

0 

4b3- 

with  ♦y  .  ^ 

-bj^  3b2 

-s 

4b, 

4b3 

0 

-bi  -by 

3b3 

0 

4b2 

4b, 

(3  X  6) 

3ai  -82 

3 

432 

0 

4a3 

and  ty  =-|j 

-ay  3ay 

■^3 

0 

-ay  -ay 

333 

0 

432 

{3  X  6) 

- 

then  B 


0  ♦y 

♦v  ♦x 

(9  X  12) 


So 


k  = 


T  T 

0  (p 

y  y 


0  4) 


T  ^  T 
y  X 


CiiQ  0^2^ 
C12Q  C22Q 


S3Q 


which  may  be  expanded  to 


1 

4  '^C,,Q(})  +4  '^C,,Q4 

^x  11  ^x  ^y  33  ^y 

1  4  '^C,  .Q4 
'  ^x  12  y 

+  4  '^QC,,4 

y  ^  33^x 

(6  X  6) 

1 

j  (6  X 

6) 

4  '^C--Q4  +  <P  '^C,,Q4 

^y  22  ^x  x  33  y 

1  T 

1  ^^22% 

+  4/33Q4j^ 

(6  X  6) 

1  (6  X 

6) 

When  k  is  partitioned  as  above,  a  small  increase  in  efficiency 
of  computing  results  as 

1/  Multiplication  by  null  matrices  is  avoided 
2/  Symmetry  of  the  stiffness  matrix  can  be  used 
3/  Storage  for  a  full  B  matrix  is  not  required. 


With 
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Tn  the  Totra-core  incxJel,  if  P-0  is  the  fiber  direction,  then 
the  node  at  the  mid-side  of  P-Q  is  never  loaded  or  supported 
out-of -piano . 


R 


Upon  generating  the  structural  stiffness  matrix,  this  node 
would,  in  general,  create  linear  dependence  in  the  stiffness 
matrix.  The  linear  dependence  could  be  removed  by  giving 
the  triangle  bending  stiffness,  but  this  would  greatly  reduce 
program  efficiency.  As  the  node  is  never  loaded,  the  simplest 
and  most  direct  solution  is  to  remove  the  node  at  the  elemental 
stage . 

This  is  done  as  follows:  The  node  to  be  removed  corresponds 
to  node  6.  So,  when  k  is  rearranged  so  that  it  corresponds 
to  the  usual  order,  the  rows  and  columns  are  rearranged  so 
that  those  corresponaing  to  1)5,  V5  occur  at  the  bottom  and 
right-hand  side  of  the  n.atrix. 
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1  •  6  •  /  {P } 


fKiil  ■  1^22!  tK2j^J 


and  we  have  the  stiffness  of  the  element  expressed  in  terms 
of  the  five  remaining  nodes  as 


[Kul 


[Kjjl  IXjjJfKjjl 


In  calculating  the  stress  matrix,  as  in  Reference  1,  we  define 
"area  coordinates:"  ■» 


Coordinates  Ci  uniquely  determine  location  P,  and  the  stress 
matrix  is 


In  the  program  this  has  been  broken  down  into  a  more  efficient 
form  as 
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where 


1X6  f,' 

4(C3b2+C2b3^4(C^b3+C3b^)> 

=  2a  ^  (4^2“!)  »  {4^2“!)  ^2 '  (4^3-l)a3,  4  (^2S2'*^^1^2^  ' 

1X6 

4(5332+^2^3),  4 {52a3+?3a2^^ 

Multiplication  by  the  C  matrix  is  carried  out  after  element 
strains  have  been  determined  in  subroutine  DEFL.  The  sixth 
node  is  reduced  out  of  the  stress  matrix  in  the  same  manner 
as  for  the  stiffness  matrix. 

(S  Tpgjj]  =  [S  T]  -  (S  T22]'^  f^22^~^^^21^ 

The  stress  matrix  is  determined  for  four  locations  in  the 
linear  strain  triangle;  Nodes  P,  Q,  and  R  and  the  triangle 
centroid. 

Quadrilateral  Element 


The  linear  strain  quadrilateral  element  used  in  the  truncated 
Tetra-core  model  is  formed  by  adding  two  linear  strain  tri¬ 
angles  from  which  the  sixth  nodes  have  been  removed. 


The  mid- side  node  between  R  and  Q  is  then  reduced  out  as  the 
top  and  bottom  mid-side  nodes  had  been,  since  it  is  never  loaded 
or  attached  to  another  element.  This  results  in  a  six-node 
linear  strain  quadrilateral.  The  stress  matrices  are  similarly 
formed,  with  stress  now  being  computed  at  the  four  corner  nodes. 

Constant  Strain  Triangle 

The  constant  strain  triangle  used  for  face  sheets  is  generated 
using  the  method  of  Przemieniecki^ .  The  local  in-plane 
coordinate  system  set  up  for  the  linear  straing  triangle 
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with  ai  =  X32f  etc.,  is  used.  A  [B]  matrix  relating  strains 
within  the  element  to  deflections  of  the  node  points  is 
defined. 


The  stiffness  matrix  [K]  is  then 
[K]  =  [C]  [B]  dV 

V 

[K]  =  At  [B]^  [C]  [B] 


Consistent  Surface  Loads  Vector 

The  uniform  in-plare  loading  acting  on  an  edge  of  the  linear 
strain  triangle  must  be  resolved  into  stress  resultants  acting 
at  the  two  corner  nodes  and  the  midpoint.  The  derivation  shown 
here  is  from  a  document  by  Tocher.^  The  total  force  on  the 
plate  is 
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The  equilibrium  equations  can  be  written  using  virtual  work 
as 


<6U,  6U.  5U-> 
14  2 


1 

/6U(0  P/3d; 


2X 

The  virtual  deflection  6U{f;)  can  be  obtained  as  (Ref.  1) 


6U  =  <5U,  61J.  6U  > 
14  2 


Thus 
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and  the  load  distribution  is  simply  two-thirds  ot  the  load 
through  the  middle  node  and  one-sixth  through  each  end  node. 


Calculation  of  Deflections 


A  solution  of  the  equation  kA  =  P  is  used  to  calculate 
deflections  (A)  from  the  stiffness  matrix  (k)  and  the  loads 
vector  (P).  A  banded  form  of  the  stiffness  matrix  is  assumed, 
in  which  only  the  nonzero  portion  of  the  stiffness  matrix, 
which  is  .’nside  the  bandwidth,  is  stored  (Figure  15).  The 
stiffness  matrix  is  stored  in  blocks  on  the  disc,  with  each 
block  being  three  rows.  Thus,  only  one  block  is  in  core  at 
any  one  time  during  generation  of  the  stiffness  matrix.  Each 
elemental  stiffness  matrix  is  generated  and  added  into  the 
overall  stiffness  matrix  using  the  direct  stiffness  method. 

The  location  of  the  block  which  each  member  of  the  elemental 
stiffness  matrix  is  added  into  is  determined  from  the  nodal 
connectivity  of  that  element.  Each  element  will  add  its 
stiffness  to  several  blocks,  thus  requiring  the  ability  to 
turn  the  needed  block  by  means  of  a  direct-access  disc. 
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Once  the  stiffness  matrix  has  been  formed,  it  is  reduced 
to  a  triangular  form  using  the  Choleski  methods  3. 

This  results  in  all  members  of  the  matrix  below  the  diagonal 


being  zero. 

i-1 

1/2 

F..  =  k.. 
11  11 

-  Z 

r=l 

j-1 

-2 

kir 

i 

=  l,n 

kij  =  kij 

k . .  =  0 

13 

-  Z 
r=l 

^ir^r  j  j  j 

i 

i 

<  j 

>  j 

where  n  is  the  order  of  the  stiffness  matrix.  The  loads 
vector  is  modified  by  the  same  factors  that  were  used  to 
triangularize  the  stiffness  matrix.  With  the  triangularized 
form  of  the  stiffness  matrix,  the  deflections  of  the  last 
node  can  be  calculated 


hnM  ’'2n^2  +  “ ^^nn^n  =  ^n 
where  JT  is  the  triangularized  stiffness  matrix. 

since  =  0 


Using  this  method,  A can  be  solved,  in  turn, 

by  back  substituting  the  previously  found  deflections  into  the 
calculations  for  the  next  higher  row.  For  example.  A,  _  can  be 
calculated 


F  .  .A  ,  +  k  ,A  =  p  . 
n-l,n-l  n-1  n,n-l  n  n-1 

where  had  been  calculated  in  the  step  before  and 

,  and  P  .  are  known 
n-1, n-1  n,n-l  n-1 


This  method  of  calculating  deflections  has  proved  to  be  much 
faster  than  the  traditional  method  of  inverting  the  stiffness 
matrix  and  multiplying  k-lp  to  get  A .  The  method  allows  a 
large  problem  to  be  solved  on  a  computer  with  small  core  space. 
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since  only  three  rows  and  three  columns  of  the  stiffness  matrix 
are  in  core  at  one  time  during  solution.  However,  it  requires 
many  accesses  to  the  disc  during  a  solution,  since  rows  and 
columns  must  be  shuffled  in  and  out  of  core  many  times.  On 
the  360/44  at  Fort  Eustis,  the  disc  space  available  is  the 
limiting  factor  on  the  size  of  problem  that  can  be  run. 

Fixities  are  applied  to  the  stiffness  matrix  by  multiplying 
the  diagonal  term  of  the  row  to  be  fixed  by  a  factor  of  10^2 , 
This  has  the  effect  of  giving  that  freedom  a  much  greater 
stiffness  than  other  freedoms,  effectively  fixing  it.  This  is 
equivalent  to  putting  a  one  on  the  diagonal  and  setting  elements 
of  that  row  and  column  to  zero,  but  it  is  more  efficient. 


OPTIMIZATION 


A  modified  form  of  the  steepest-descent  opt-_mi?,ation  method 
is  used  to  find  a  minimum  weiqht  Tetra-core  design  for  a  given 
set  of  loads^.  Variables  that  can  be  optimized  are 
tetrahedron  geometry  (side  length,  height,  theoretical  height, 

Yl.  offset,  Y2  offset,  THETA)  and  the  thicknesses  of  each  leg. 

See  Figures  4  and  5  for  a  description  of  the  tetrahedron  geo¬ 
metry.  Optimdzation  controls  have  been  set  up  so  that  any 
combination  of  variables  can  be  held  constant  during  optimi¬ 
zation.  Also,  any  lihree  variables  can  be  forced  to  be  the 
same  during  optimization;  i.e..  Vertical  and  Skew  A  and  B  legs 
can  be  kept  to  the  same  thickness  during  optimization.  Input 
Card  Type  2  is  used  to  specify  which  vari  ihles  v.'ill  be  optimized 

The  optimization  method  will  be  illustrated  using  a  two- 
dimensional  example  so  that  the  variation  of  ch  sign  variables 
can  bo  plotted  graphically.  For  tlie  exiimple,  assume  that  a 
Totra-core  run  is  being  made  in  which  all  variables  are  held 
constant  except  for  the  vertical  leg  thickness  and  the  skew 
leg  thickness,  and  that  the  Skew  A  and  B  legs  will  have  the 
same  thicknesses.  Using  this  model,  a  design  space  can  be 
plotted  for  a  given  set  of  loads  in  wiiich  the  X  and  Y  axes 
correspond  to  the  variation  in  vertical  and  skew  leg  thicknesses 
as  shown  in  Figure  16.  For  any  point  in  the  design  space,  the 
Tetra-core  element  will  have  a  given  weight  and  margin  of  safety 
A  series  of  curve.s  of  constant  v;eighl  and  another  series  of 
curves  of  constant  margin  of  safety  can  be  dravm  as  shown  in  the 
figure.  The  margin-of- safety  curve  for  a  margin  of  safety  of 
zero  is  used  to  divide  the  design  space  into  a  feasible  region, 
in  v;hich  the  margin  of  safety  is  greater  than  .:ero,  and  an 
unfeasible  region,  with  a  margin  of  safety  less  than  zero.  The 
optimization  program  does  not  know  the  location  cr  shape  of 
these  curves,  except  for  the  location  of  the  present  design,  but 
it  can  calculate  the  gradient  of  these  curves;  that  is,  the 
direction  of  travel  in  desion  space  which  will  result  in  the 
maximum  change  in  either  the  weight  or  the  m.argin  of  safety. 

The  gradient  can  be  calculated  by  changing  each  variable, by  a 
small  amount,  while  holding  all  ether  variables  constant  and 
calculating  the  weight  and  margin  of  safety  for  the  new  design 
point.  This  is  done  separately  for  all  variables,  resulting  in 
a  set  of  partial  derivatives: 
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where  and  refer  to  the  design  point  at 

the  start  of  calculation  of  derivatives,  before  changing  any 
thickness.  Then  these  sets  of  partial  derivatives  can  be  used 
to  calculate  the  required  gradients. 
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M  =  number  of  constraints 
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These  gradients  are  actually  the  direction  cosines  of  the 
normal  to  each  curve  at  that  design  point. 

For  example,  assume  that  the  starting  design  point  is  chosen  at 
Point  A  in  Figure  17  and  that  the  weight  and  margin  of  safety 
have  been  calculated  for  this  point.  As  yet,  no  partial 
derivatives  have  been  calculated,  so  the  thickness  will  be 
decreased  in  steps  in  the  same  ratio  to  get  t.o  the  constraint 
curve  (MS=0) .  In  the  general  case,  the  program  changes  only  the 
thicknesses  in  first  reaching  the  constraint  curvi’,  leaving  t!;e 
side,  height,  etc.,  dimensions  constant.  Then  the  program 
calculates  the  gradients  of  the  weight  and  the  margin-of-saf ety 
curves  as  previously  described.  There  is  a  different  margin  - 
of-safety  curve  for  each  leg,  since  changing  the  thickness  of 
one  leg  will  change  the  margin  of  safety  in  that  log  more  than 
in  another  leg.  This  requires  the  calculation  of  a  gradient  of 
the  margin  of  safety  in  each  leg,  as  shown  in  Figure  IB.  A 
composite  constraint  gradient  is  computed  using  the  gradient 
to  each  constraint  curve  and  the  margin  of  safety  in  each  leg. 

{y}  =  (1-MS^)^^  +  {B2)  (1-MS2)^^  +...+  (1-MSj^)^^ 

The  proportion  in  which  each  constraint  gradient  is  added  into 
the  composite  constraint  gradient  is  determined  by  the  distance 
of  the  design  point  from  that  constraint. 


In  Figure  18  for  case  A,  in  which  the  design  point  is  at  the 
intersection  of  two  constraints,  each  gradient  is  added  in 
equally,  resulting  in  a  new  direction.  In  case  B,  in  which 
the  design  point  is  on  one  constraint  and  far  from  the  other, 
the  gradient  of  the  constraint  which  the  design  is  on  will 
have  the  major  influence  on  the  new  direction.  This  composite 
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Figure  17.  TwoOimeiuional  Optimization  Example. 
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Figure  18.  Composite  Constraint  Gradient. 
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constraint  gradient  is  used  by  thc’  program  to  make  a  step  in 
the  direction  of  maximum  change  in  margin  of  safety,  to  Point  C 
in  Figure  17.  The  step  is  caicalated: 


where  X  is  the  step  length.  Next,  the  design  variables  are 
decreased  along  the  gradient  to  the  constant  weight  curve  to 
react,  the  starting  weight  (Point  D)  .  This  cycle  B-C-D  results 
in  a  sidestep  av/ay  from  the  constraint  curve  in  the  direction 
of  maximum  increase  in  margin  of  safety,  while  maintaining 
a  constant  weight.  The  margin  or  safety  is  calculated  at 
Point  D.  If  the  margin  of  safety  is  now  greater  than  it  was 
at  Point  B,  then  the  bucket  in  the  curve  has  not  been  reached 
and  another  sidestep  is  taken  (D-F) .  Jtirgin  of  safety  is 
recalculated  and  another  sidestep  is  taken  (F-E) .  When  the  margin 
of  safety  is  calculated  at  Point  K,  it  is  found  to  be  less  than 
at  Point  F.  Therefore,  the  bucket  in  the  constraint  curve  has 
been  passed  and  a  backward  step  is  taken  (E-F) .  Now  the  vari¬ 
ables  are  decreased  along  the  gradient  to  she  constant  weight 
curve  (this  is  known  as  the  direction  of  steepest  descent)  in 
steps, to  reach  the  constraint  curve  (Point  G) .  This  scries  of 
operations  (B-G)  is  called  one  cycle  in  the  program.  At  the 
completion  of  one  cycle  (Point  G) ,  gradients  are  recalculated, 
the  step  length  is  reduced,  and  another  series  of  sidesteps 
is  begun.  The  number  of  optimization  cycles  to  be  run  is  input 
as  "NSTP"  on  card  type  one.  If  a  zero  or  blank  is  input  for 
"NSTP,"  the  program  v/ill  set  "NSTP"  -  3.  If  the  program  makes 
ten  sidesteps  (B-D,  D-F,  etc.)  without  reaching  the  bucket, 
the  program  will  descend  to  the  constraint  curve,  recalculate 
gradients,  and  begin  sidestepping  again  with  a  larger  step 
length  to  reach  the  optimum  in  a  fewer  number  of  steos.  At 
the  end  of  each  cycle  (B-G),  the  v;C'ight  of  the  design  is  checked. 
If  the  new  weight  is  less  than  the  weight  at  the  end  of  the 
previous  cycle,  the  variables  are  stored  as  a  new  optimxun  design. 
Generally,  the  program  v/ill  find  a  series  of  designs,  each  lighter 
than  the  previous  one.  However,  it  is  possible  that  the  optimum 
will  be  found  on  the  first  or  second  cycle,  and  the  rest  of  the 
cycles  v.'ill  result  in  heavier  designs.  This  makes  it  necessary 
to  store  the  lightest  design  found,  since  the  program  may  not 
return  to  it. 

In  general,  the  constraint  and  constant  weight  curves  are  not 
two  dimensional  but  are  N-dimensional  hypersurfaces,  N  being 
the  number  of  variables.  Figure  19  illustrcates  the  set  of 
constraint  surfaces  for  a  three-dimensional  problem.  Generally, 
several  optimum  designs  will  be  possible  for  a  given  set  of 
variables,  and  the  program  will  only  find  the  one  that  is 
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nearest  to  the  input  starting  design.  In  Figure  20,  if 
optimization  is  begun  from  Point  A,  an  optimum  design  with 
a  weight  of  0.9  pound/square  inch  will  be  found,  while  starting 
at  Point  B  will  result  in  a  weight  of  0.08  pound/square  inch. 

Thus  it  is  necessary  to  make  several  optimization  runs  to 
determine  if  the  "global"  or  most  optimum  design  has  been  found. 

The  steepest  descent/sidestep  optimization  method  can  be  thought 
of  as  a  way  of  getting  down  a  mountain  pass  to  the  lowest  alti¬ 
tude  possible.  The  sides  of  the  mountain  represent  constraints, 
and  lines  of  constant  altitude  represent  lines  of  constant  weight. 
Once  a  given  altitude  has  been  reached  by  the  design  point,  the 
altitude  cannot  be  increased — only  decreased. 

A  plot  of  the  weight  versus  number  of  cycles  for  a  typical 
case  is  shown  in  Figure  21.  For  this  run  a  true  Tetra-core 
flat  plate  was  optimized  with  the  side  length  held  constant, 
the  height  and  vertical  leg  thickness  allowed  to  vary,  and  the 
Skew  A  and  B  leg  thicknesses  held  the  same.  The  largest  change 
in  weight  was  in  the  first  cycle,  indicating  that  the  optimum 
region  was  found  quickly,  with  the  rest  of  the  cycles  being 
taken  up  with  searching  the  immediate  region  to  find  the  exact 
location  of  the  optimum. 


NONLINEAR  ANALYSIS 


An  analysis  of  the  effect  of  the  nonlinear  stress-strain 
behavior  of  the  fiber-glass-epoxy  composite  used  in  Tetra-core 
legs  can  be  conducted  with  this  program.  Stress-strain  curves 


of  the  composite,  determined  by  cutting  out  legs  from  the  panel 
and  testing  them  separately,  are  used  as  input, as  is  the  loading 
increment  to  be  applied^.  First, the  program  calculates 
the  stresses  and  strains  in  each  plate,  using  the  tangent  moduli 
and  Poisson's  ratios  input  on  Card  Type  4  and  the  loads  input 
on  Card  Types  12  and  13.  -11160,3  new  modulus  is  calculated  for 

each  plate, using  the  stress -strain  curves  input  for  each  leg 


and  the  strains  already  calculated.  As  an  example,  suppose 
that  the  strain  in  the  fiber  direction  (e^)  of  one  plate  had 
been  +  O.OOli  at  the  end  of  the  first  load  increment,  and  that  the 
tension  stress-strain  curve  shown  in  Figure  22  had  been  input  on 
Card  Typo  8  as  a  series  of  points  for  the  leg  which  this  plate 
was  in,  then  the  program  takes  the  strain  of  0.0011  and  usas  a 


linear  curve  fitting  method  to  determine  the  stress  corresponding 
to  the  given  strain.  This  stress  is  used  to  compute  a  secant 
modulus  tc  be  used  by  the  plate  for  the  second  iteration: 


sec  1  1 

The  secant  modulus  is  determined  in  the  longitudinal,  transverse, 
and  shear  directions  from  the  stress-strain  curves  input  for 
those  directions.  A  different  stress-strain  curve  is  input  for 
tension  and  compression  for  the  longitudinal  and  transverse 
directions, and  the  correct  one  is  chosen  by  referring  to  the  sign  of 
the  strain  in  that  direction. 

Once  moduli  and  Poisson's  ratio  in  each  plate  have  been  determined, 
the  loads  are  increased  by  the  amount  input  on  Card  Types  12 
and  13  and  a  new  stress  analysis  is  run.  Then  the  stresses  and 
strains  in  each  plate  are  used  to  determine  a  new  set  of  secant 
moduli  and  Poisson's  ratio  for  each  plate.  The  loads  are  increased 
again  and  the  stresses  are  recalculated.  This  loading  in  increments 
continues  until  the  number  of  steps  specified  on  Card  Type  1  as 
"NSTP"  are  completed.  As  the  loads  are  increased,  some  plates  will 
reach  their  failing  strain  in  some  direction.  When  this  happens, 
the  secant  modulus  in  that  direction  is  reduced  to  a  small  number, 
which  unloads  that  plate  without  allowing  the  nodes  to  which  it 
is  attached  to  become  free,  as  would  happen  if  the  modulus  were 
reduced  to  zero.  On  the  next  loading  increment,  strains  in  the 
adjoining  plates  will  be  higher  than  would  be  caused  by  the  increase 
in  end  load  only,  since  the  loads  in  the  failed  plate  have  been 
dumped  into  them.  This  will  result  in  the  failure  of  more  plates  in 
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Figure  22. 


Method  Used  To  Determine  Moduli  for  Plastic  Analysis 


the  same  area.  Finally,  enough  plates  will  have  failed  that 
there  is  no  load  path  through  the  structure,  resulting  in 
complete  failure.  Figure  23  illustrates  a  nonlinear  analysis 
run  on  a  Tetra-core  panel  with  several  plates  removed  to 
simulate  the  effect  of  a  hole.  Initially,  all  plates  except 
those  that  have  been  removed  are  intact,  but  as  the  load 
is  increased, several  plates  around  the  hole  fail;  and  finally 
at  a  load  of  3,000  pounds/inch,  total  failure  occurs,  if  a 
lower  increment  of  load  had  been  used,  the  spread  of  failure 
around  the  hole  would  have  been  more  visible. 
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Figure  23  Example  of  Nontincar  Anaiv!,i$, 
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BUCKLING  ANALYSIS 

A  buckling  analysis  is  performed  on  the  plate  elements  in 

each  leg  making  up  the  Tetra-core  panels  in  the  subroutine 

BUKL.  No  equation  was  available  to  calculate  the  compression 

buckling  load  of  an  orthotropic  triangular  plate  with  one  edge 

free,  loaded  in  compression.  Therefore,  it  was  decided  to 

use  the  gompression  buckling  equation  (developed  by  Lackman 

and  Ault  )  for  a  rectangular  orthotropic  plate  with 

one  edge  free  and  three  edges  simply  supported.  The  length 

of  the  free  edge  on  the  triangle  was  used  as  the  length  of 

the  rectangle,  with  the  tetrahedron  height  used  as  a  rectangle 

width.  This  simplification  is  expected  to  result  in  conservative 

values  for  the  buckling  load,  due  to  the  greater  restraint  imposed 

on  the  panel  by  the  sides  of  the  triangle  compared  with  the  sides 

of  a  rectangle.  The  equation  used  to  calculate  the  buckling  load 

is 
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a  =  tetrahedron  side  length 
b  =  tetrahedron  height 
t  =  leg  thickness 

An  iteration  method  is  used  to  find  the  lowest  value  of  Nx 
which  will  satisfy  the  equality.  This  calculation  is  made 
for  the  Vertical,  Skew  A,  and  Skew  B  leg.s.  If  face  sheets  are 
applied  to  the  Tetra-core  panel,  then  all  edges  of  the  triangle 
will  be  simply  supported  and  the  equation  of  a  simply  supported 
rectangle  is  used'. 

2  1/2 

N^-2(I)  ((DjD^)  ♦  D3I 


This  equation  is  used  for  vertical.  Skew  A  and  B  legs,  and 
upper  and  lower  face  sheets  when  face  sheets  are  applied.  If 
a  more  accurate  buckling  equation  becomes  available,  it  can 
easily  be  inserted  into  the  buckling  subroutine'. 


MARGIN  OF  SAFETY 


The  margin  of  safety  is  calculated  for  each  plate  element  using 
the  von  Mises-Hill-Tsai  equation^. 
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01,  02,  and  T12  are  the  stresses  as  in  the  local  plate  element 
coordinate  system,  with  ai  always  in  the  fiber  direction.  Fx, 

Fy,  and  Fxy  are  the  allowable  stresses,  with  Fx  being  the 
allowable  stress  in  the  fiber  direction,  Fy  the  transverse 
allowable,  and  Fxy  the  shear  allowable.  Tension  or  compression 
allowable  stresses  (as  input  on  card  type  10)  are  used,  depending 
on  the  sign  of  the  stress.  A  margin  of  safety  is  calculated  for 
the  four  locations  in  the  plate  for  which  stresses  have  been 
calculated,  with  the  lowest  margin  being  printed  out  forthat 
plate.  The  margin  of  safety  in  buckling  is  also  calculated, 
using  the  buckling  stress  calculated  for  the  leg  in  which  this 
plate  is  found. 


The  buckling  margin  is  calculated  only  if  one  of  the  a]|s  is 
negative,  indicating  a  compression  load  in  the  plate.  The 
lowest  margin  of  safety  found  in  each  leg  for  that  load  case  is 
printed  after  printing  of  the  plate  element  stresses. 
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CALCULATION  OF  ULTIMATE  LOADS 


After  the  margin  of  safety  in  each  plate  has  been  determined, 
the  minimum  margin  of  safety  found  in  any  plate  is  used  to 
extrapolate  the  given  applied  load  to  determine  the  load  at 
which  the  model  would  have  failed. 
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where  Y^s  is  the  minimum  margin  of  safety  found  in  any  plate. 
J^load/  Yloadf  etc.,  are  the  input  loads  from  card  types  12  and 
13.  The  power  of  0.5  may  have  to  be  changed  to  agree  with 
tests  when  they  become  available. 

CALCULATION  OF  STIFFNESSES 


An  equivalent  modulus  in  the  x  direction  is  calculated  whenever 
a  uniform  load  is  applied  in  that  direction.  This  is  done  by 
taking  two  x  deflections  on  opposite  ends  of  the  panel  and 
adding  them  to  determine  the  total  change  in  length  of  the  panel 
(AX) .  This  AX  is  divided  by  the  original  panel  length  to  get 
an  equivalent  strain  The  uniform  load  on  the  x  face  of 

the  panel,  divided  by  the  height  of  the  face,  gives  an  equivalent 
stress  Ox„,s.  The  equivalent  modulus  in  the  x  direction  is  then 


E  =  — ^ 
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An  equivalent  modulus  in  the  y  direction  is  calculated  in 
the  same  way  whenever  a  uniform  load  is  applied  to  the  y  face. 

OUTPUT  DESCRIPTION 

Primary  output  of  the  Tetra-core  analysis  program  consists  of 
the  finite-element  model  generated  by  the  data  generator  over¬ 
lay  and  the  deflections  and  stresses  calculated  by  the  stress 
analysis  overlay.  If  the  optimization  option  is  used,  additional 
printout  is  generated  during  each  step  of  optimization. 

The  printed  output  from  a  standard  stress  analysis  run  starts 
with  the  input  data  (number  of  load  cases,  size  of  tetrahedron, 
leg  thicknesses,  etc.);  the  printout  of  calculated  data  starts 
with  the  number  of  node  points  and  pJate  elements  generated 
by  the  data  generator.  This  node  point  number  does  not  include 
the  nodes  at  the  midpoint  of  the  triangle  or  quadrilateral 
which  are  added  in  the  stress  analysis  overlay.  The  number  of 
load  cases  is  also* printed.  If  a  cylinder  is  generated,  the 
radius  to  the  outer  surface  is  printed.  If  an  airfoil  is 
generated,  the  chord  length  is  printed.  The  X,  Y,  Z  coordinates 
of  the  generated  nodal  points  are  printed.  Nodal  loads  are 
printed  for  each  load  case  along  with  the  original  node  number 
and  the  new  node  numbers  that  include  mid-side  nodes  added  by 
each  linear  strain  triangle  and  quadrilateral.  Each  linear 
strain  element  used  in  the  model  adds  a  new  node  point  at  the 
midpoint  of  two  sides  of  the  element  where  it  is  attached  to 
another  element,  which  must  be  accounted  for  by  increasing  the 
size  of  the  load  vector  and  the  stiffness  matrix.  The  new 
node  numbers  are  generated  from  the  node  numbers  of  the  vertices 
on  each  side  of  the  new  node  (Figure  24).  The  new  node  number 
is  generated  by  multiplying  the  larger  of  the  two  nodes  by 

I, 000  and  adding  to  the  smaller,  resulting  in  a  unique  new  node 
number.  For  example,  if  a  new  node  is  to  be  generated  between 
nodes  5  and  3, as  shown  in  the  figure,  the  resulting  number 

would  be  5003.  This  node  is  then  added  into  the  stiffness  matrix 
and  load  vector  behind  the  location  of  the  larger  of  the  two 
vertex  nodes.  The  added  mid-side  nodes  for  a  typical  case  are 
shown  in  Figure  26. 

Plate  properties  are  printed.  The  node  numbering  sequence 
for  each  plate  in  the  original  node  point  order  and  as 
renumbered  is  printed.  Nodes  listed  under  the  headings  I, 

J,  K,  and  L  give  the  original  plate  connectivity.  Nodes  listed 
under  LP,  LQ,  LR,  LS,  LT,  and  LU  give  the  connectivity  using  the 
new  nodal  numbering  sequence.  LS  and  LT  represent  new  mid¬ 
side  nodes  for  the  triangle.  Plate  moduli  and  plate  thickness 
are  printed  along  with  the  node  numberings.  EX,  EY,  and  EXY 
have  been  divided  by  the  factor  l-Uxy'Myx  'Jse  in  the  element 
stiffness  generator,  EX  is  the  Young's  modulus  in  the  direction 
of  nodes  I  and  J,  and  EY  is  at  an  angle  of  90  degrees  to  EX. 
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LINEAR  STRAIN  TRIANGLE 


Out  of  Element  Stiffness  Matrix 


Y  LINEAR  STRAIN QUADRIUTERAL 


Been  Reduced  Out 


CONSTANT  STRAIN  TRIANGLE 


Local  )^Axis  Is  in  Direction  of  Nodes  l-^J 

Local  Y-Axis  Is  Perpendicular  to  X-Axis  in  Direction  of  Node  K 


X 


Figura  24  Etamants  Uiad  in  Tatra-Cora  Program . 
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Numbar  of  Vertex  Node*  °  72 
Number  of  Mid-Side  Nodes  =  96 

Total  Number  of  Nodes  =  1 68 


.u 


•  Original  Vertex  Node 
Added  Mid-Side  Node 


♦  Y 


Figure  26i  Added  Mid-Side  Nodei,  True  Tetra-Core  Flat  Plate 
*  6,  Ly  *  6* 
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The  1/2  bandwidth,  in  nodes,  required  to  store  the  stiffness 
matrix  is  printed  after  the  elements  have  been  renumbered 
and  printed.  The  storage  space  required  for  the  overall 
stiffness  matrix  can  be  calculated  by  multiplying  the  number 
of  nodes  times  3  times  the  1/2  bandwidth  times  3.  If  this 
number  added  to  the  amount  of  space  required  to  store  the 
loads  vectors  is  greater  than  the  available  disc  space,  then 
the  problem  will  have  to  be  reduced  in  size. 

Nodal  deflections  and  plate  stresses  are  printed  for  each  load 
case  separately.  Deflections  are  printed  for  each  node  with 
the  node  points  numbered  in  the  original  sequence  (1,  2,  3, 

3001,  4,  4002,  etc.).  Stresses  are  printed  for  four  locations 
in  the  linear  strain  element  used  in  Vertical  and  Skew  A  and 
Skew  B  legs,  and  for  the  triangle  center  for  the  constant  strain 
triangle  used  in  face  sheets.  The  connectivity  of  each  triangle 
is  printed  in  the  original  numbering  sequence.  For  a  linear 
strain  triangle,  the  stresses  are  printed  at  node  I,  node  J, 
node  K,  and  at  the  center  of  the  triangle,  in  that  order.  For 
a  linear  strain  quadrilateral,  which  is  used  only  in  the 
truncated  Tetra-core  model,  stresses  are  printed  at  nodes  I,  J, 

K  and  L.  The  margin  of  safety,  calculated  using  the  Hill-Tsai 
failure  criterion,  in  each  element  and  the  leg  which  that  element 
is  in  (Vertical  leg  =  leg  1,  Skew  A  =  leg  2,  Skew  B  =  leg  3, 
lower  face  =  leg  4,  upper  face  =  leg  5)  are  printed  also.  The 
buckling  stress  of  a  triangle  in  each  leg  is  printed. 

After  all  the  plate  stresses  in  a  load  case  have  been  printed, 
the  minimum  margin  of  safety  found  in  each  leg  is  printed. 

Then  the  minimum  margin  of  safety  is  used  to  extrapolate  the 
original  loads  to  predict  a  load  case  with  the  same  ratio 
between  each  load  that  would  have  a  margin  of  safety  of  zero. 

If  the  model  has  been  loaded  in  the  X  or  Y  direction,  the 
deflections  of  two  node  points  on  opposite  ends  of  the  model 
are  used  to  calculate  an  equivalent  modulus.  Then  printing 
of  the  deflections  and  stresses  for  the  next  load  case  starts. 

If  the  nonlinear  option  is  specified,  a  series  of  load  increments 
will  be  applied  until  overall  failure  occurs.  Nodal  coordinates 
and  load  vectors  will  be  printed  out  only  for  the  first  incre¬ 
ment,  since  the  coordinates  do  not  change  and  the  succeeding  load 
vectors  will  be  integer  multiples  of  the  first  load  vector.  The 
load  intensity  (i.e.,  load  in  X  direction  in  pounds/inch)  is 
printed  at  the  beginning  of  each  load  increment.  Plate  element 
moduli  are  printed  at  the  start  of  each  load  increment  to  allow 
an  inspection  of  the  moduli  in  each  plate  to  determine  the  extent 
of  degradation  in  the  panel.  As  each  plate  fails  in  a  given 
direction  (longitudinal,  transverse,  or  shear),  the  modulus  in 
that  direction  is  reduced  to  near  zero.  Final  failure  occurs 
when  all  plates  have  failed.  Deflections  and  stresses  are 
printed  out  for  each  load  increment. 
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If  the  hole  option  is  specified^  the  moduli  each  triangle 
that  connects  to  a  node  in  the  center  of  thi>  plate,  cylinder, 
or  airfoil  are  set  to  zero.  This  node  is  normally  determineJ 
by  the  program  but  can  be  input  if  desired.  The  identity  of 
the  affected  node  is  printed  out  in  either  case. 

If  the  optimization  option  is  specified,  the  coordinates,  loads, 
elements,  deflections,  and  stresses:  are  printed  out  for  the 
input  configuration  and  for  the  final  optimized  configuration 
only.  Intermediate  output  during  each  optimization  step 
consists  of  the  current  design  configuration  and  the  margin  of 
safety  in  each  leg.  When  partial  derivatives  are  being  calcu¬ 
lated,  the  variable  that  is  being  changed  is  printed  out  as 
"IK".  Direction  cosines  of  the  constant  weight  surface, 
direction  cosines  of  the  constraint  surface  of  each  leg,  margin 
of  safety  in  each  leg,  and  direction  cosines  of  the  composite 
constraint  surface  are  printed  at  the  start  of  each  cycle  of 
sidestepping.  Weight  of  the  design  is  checked  at  the  end  of 
each  sidestepping  cycle  to  determine  if  a  new  optimum  has  been 
reached.  If  so,  the  design  variables  are  printed  and  optimi¬ 
zation  continues.  When  the  input  number  of  sidestep  cycles 
"NSTP"  have  been  completed,  the  optimization  process  is  stopped 
and  the  deflections  and  stresses  for  the  lowest  weight  design 
that  was  found  are  printed. 


64 


SAMPLE  PROBLEMS 


Examples  of  input  for  several  types  of  analyses  will  be  given 
for  reference.  For  a  nonlinear  analysis,  the  set  of  cards  shown 
in  Figure  26  will  be  required.  Three  sets  of  cards  of  Types  8 
and  9  are  required  to  input  stress-strain  curves  for  each  leg. 
The  input  cards  required  for  an  airfoil  analysis  are  shown  in 
Figure  27.  An  example  of  the  input  for  an  optimization  run  is 
shown  in  Figure  28. 

The  output  from  an  optimization  run  is  shown  in  Appendix  II. 
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Figure  2&  Input  Data  for  Nonlinear  Analysis. 


Figure  27.  Input  Data  for  Airfoil  Analysis. 


Card  Type  1 


Figur*  28.  Input  Data  for  Optimization  Run. 
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COMPARISON  WITH  UNIVERSITY  OF  KANSAS  THSTS 


Computer  runs  have  been  made  to  compare  to  several  of  the 
tests  run  at  the  University  of  Kansas^.  The  results  of  this 
comparison  are  not  very  meaningful,  however,  because  the 
basic  properties  of  the  material  making  up  the  legs  were 
not  specified  in  the  test  document.  Information  in  the  test 
document  was  limited  to  side  length  and  tetrahedron  height 
and  the  fiber  used,  which  was  S  glass.  The  program  requires 
elastic  moduli  and  thickness  for  each  leg.  It  was  assumed 
that  the  material  had  the  following  moduli  and  thickness: 

=  4,000,000  psi 
Ey  500,000  psi 
^XY  “  500,000  psi 


t  =  0.030  in. 

Using  these  as  input  with  a  side  length  of  C.6123  inch  and  a 
height  of  0.4  inch  (comparable  to  specimen  1  in  Reference  0) 
produced  an  Ex  of  281,000  psi  and  an  Ey  of  284,000  psi.  The 
predicted  initial  failure  load,  using  a  linear  stress  analysis, 
was  789  pounds.  The  cause  of  this  failure  was  a  high  transverse 
tension  stress  in  some  Skew  A  and  B  plates  which  v;ould  not  have 
resulted  in  final  failure.  If  the  criterion  of  failure  had  been 
fiber  stress,  the  predicted  ultimate  load  would  have  beer.  3,300 
pounds  for  a  fiber  stress  of  600,000  psi.  This  compares  to  an 
Ex  of  363,053  psi,  an  Ey  of  303,198  psi,  and  an  ultimate  load 
of  5,800  pounds  from  the  test.  How^ever .  a  comparison  between 
the  analysis  and  test  would  be  risky,  since  material  properties 
of  the  specimens  used  for  the  test  were  guessed  at.  The  low 
predicted  load  of  789  pounds  v/as  caused  by  failure  of  the 
Skew  A  and  B  legs  in  the  transverse  direction. 

An  analysis  run  with  a  configuration  comparable  to  that  of 
Specimen  3  (HT  =  0.5,  .side  =  0.75)  was  run,  resulting  in  an 
Ex  of  293,400  psi  that  compares  to  the  test  result  of  Ev  = 
288,150  psi.  This  appears  to  be  good  agreement,  assuming 
that  the  input  moduli  and  thickness  were  close. 
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MODEL  CONSTRAINTS 


The  Tetra-core  analysis  program  is  dimensioned  for  1100  plate 
elements,  300  vertex  nodes,  and  600  total  nodes.  A  problem  of 
this  size  cannot  be  run  on  the  IBM  360/44  due  to  a  lack 
of  disc  space.  However,  the  program  can  be  converted 
to  a  larger  360,  such  as  a  360/65,  very  easily,  should  one 
become  available.  This  can  be  done  by  removing  the  call  load 
cards  which  are  used  to  load  a  new  overlay  before  referring 
to  it.  These  cards  are  not  necessary  on  the  operating  system 
used  by  the  larger  IBM  machines.  Also,  job  control  cards 
must  be  changed.  With  a  larger  machine  it  will  be  possible 
to  increase  the  size  of  problem  that  may  be  run  to  600,  or 
greater  vertex  nodes.  This  can  be  done  by  increasing  the 
X,  Y,  Z  and  UX,  UY,  UZ  dimensions  in  the  generation  and  stress 
analysis  subroutines.  The  dimensions  of  variables  A  and  DIAG 
in  subroutine  STIFF  must  be  increased.  In  subroutine  SOLPAC, 
variables  DIAG,  AI,  AJ,  B,  and  X  must  be  increased  in  size. 

In  subroutine  DEFL, variables  IDO  and  B  must  be  increased  in  size. 

In  converting  to  a  CDC  6600  machine,  it  will  be  necessary  to 
put  in  OVERLAY  and  CALL  OVERLAY  cards  for  each  overlay.  Direct 
access  read/write  cards  used  for  the  IBM  360  will  not  work  on 
the  CDC  6600.  Direct-access  read/write  commands  will  have  to 
be  rewritten  for  the  particular  CDC  6600  on  which  the  program 
is  to  be  run,  since  they  vary  from  installation  to  installation. 

The  truncated  Tetra-core  model  can  generate  the  full  range  of 
Tetra-core  geometries,  from  a  true  Tetra-core  to  a  truncated 
Tetra-core  with  vertical  sides.  However,  both  the  true  Tetra- 
core  and  the  vertical-sided  Tetra-core  models  will  have  several 
nodes  which  are  connected  to  each  other  by  plates  generated  for 
the  same  X,  Y,  Z  coordinates.  This  will  cause  a  linear  depend¬ 
ence  in  the  stiffness  matrix,  since  these  nodes  will  essentially 
have  the  same  stiffnesses;  that  is,  their  rows  and  columns  will 
be  identical.  This  will  result  in  the  matrix  being  singular, 
and  no  solution  can  be  found.  This  problem  will  also  result 
when  the  connected  nodes  are  within  a  small  distance  of  each 
other. 

The  program  does  not  compute  nodal  reaction  force,  since  this 
would  require  more  disc  storage  space  than  is  available  on  the 
IBM  360/44  at  Fort  Eustis.  If  the  program  is  converted  to  a 
larger  machine,  reactions  can  be  calculated  by  storing  a  copy 
of  the  merged  stiffness  matrix  and  then  multiplying  it  by  the 
calculated  deflections:  P  =  KA. 
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COMPUTER  PROGRAM 


Input  Formats  for  Tetra-Core  Analysis  Program 


CARD 

TYPE  VARIABLE  FIELD  COMMENT 

1  NLOD  I  (1-4)  Number  of  load  cases  to  be  run.  If 

"lOPT"  =  1  or  "NLIN"  =  1  then  NLOD 
must  =  1. 

lOPT  I  (5-8)  If  "lOPT"  =  0,  no  optimization  will 

be  done.  If  "lOPT"  =  1,  the  xetra- 
core  element  will  be  optimized  for 
given  load  case. 


NLIN  I  (9-12'  "NLIN"  =  D,  Standard  linear 

stress  analysis  will  be  run.  If 
"NLIN"  =  1,  a  nonlinear  analysis  will 
be  run  using  the  input  load  case  as 
the  load  increment  and  "NSTP"  as  the. 
number  of  iterations. 


ITYP  I  (13-16)  If  "ITYP"  =  1,  flat  plate  analysis 

If  "ITYP"  =  2,  cylinder  analysis 
If  "ITYP"  =  3,  airfoil  analysis 

IFACE  I  (17-20)  If  "IFACE"  =  0 ,  no  face  sheets  will 

be  added  to  the  Tetra-core. 

If  "IFACE"  =  1,  face  sheets  will  be 
added. 


IHOLE  I  (21-24)  If  "IHOLE"  =  0,  a  standard  analysis 

will  be  run.  .  If  "IHOLE"  =  1,  all 
plates  connected  to  a  grid  point  in 
the  center  of  the  panel  will  be  de¬ 
fined  with  zero  stiffnesses,  giving 
the  effect  of  a  hole  in  the  panel. 

If  "IHOLE"  >  1,  all  plates  connected 
to  grid  point  "IHOLE"  will  have  zero 
stiffness,  allowing  a  hole  to  be 
inserted  where  desired. 

NSTP  I  (25-28)  "NSTP"  determines  the  number  of  itera 

tions  to  be  made  in  a  nonlinear  analy 
sis,  or  the  number  of  optimization 
cycles  to  be  made  in  an  iptimization 
run.  The  default  option  for  "NSTP" 
is  3. 
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CARD 

TYPE  VARIABLE  FIELD  COMMENT 


ITET  I  (29-32)  If  "ITET"  =  0,  standard  tetrahedrons 

will  be  generated.  If  "ITET"  =  1, 
truncated  tetrahedrons  will  be 
generated. 


ITEST  I  (33-36)  If  "ITEST"  =1,  the  element  and 

merged  stiffness  matrices  will  be 
printed  out.  If  "ITEST"  =  0 ,  no 
printout  of  matrices. 

ILOD  I  (37-40)  If  "ILOD"  =  0,  the  program  will 

automatically  generate  nodal  loads 
and  fixities.  If  "ILOD"  >  0,  nodal 
loads  and/or  fixities  are  input  on 
Card  Types  14  and  15,  as  shown  below. 


ILOD  Comment 


1 

2 

3 


Nodal  loads  are  input  on  Card  Type  14. 

Nodal  fixities  are  automatically  generated. 

Nodal  loads  are  input  on  Card  Type  14 . 

Nodal  fixities  are  input  on  Card  Type  15. 

Nodal  loads  are  automatically  generated. 
Nodal  fixities  are  input  on  Card  Type  15. 
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CARD 


TYPE  VARIABLE 

FIELD 

COMMENT 

2  ISA 

I 

(1-4) 

Card  Type  2  is  input  only  if  "OPT"  = 
1.  "ISA",  "ISB",  "ISC"  are  used  to 

ISB 

I 

(1-5) 

set  variables  equal  to  each  other 
during  optimization.  11  cards  of 
type  2  are  input. 

ISC 

I 

(9-12) 

During  optimization  the  program,  will 
normally  change  each  variable  inde- 

pendently.  ISA,  ISB,  ISC  can  be  used 
to  keep  several  variables,  such  as 
leg  thicknesses,  the  same  during  opti¬ 
mization.  ISA  can  also  be  used  to 
hold  variables  constant  during  opti¬ 
mization  by  inputting  ISA  =  12  for 
that  variable.  One  Card  Type  2  is 
input  for  each  variable,  resulting  in 
eleven  cards  of  type  2.  The  first 
Card  Type  2  is  used  to  control  the 
side  length,  the  second  controls 
height,  the  third  controls  vertical 
leg  thickness,  etc.,  as  shown  below. 

If  ISA  =  2  on  the  first  card  type  2 
and  ISA  =  1  on  the  second  card  type 
2  then  the  first  and  second  variables 
(side  length  and  height)  will  be  kept 
the  same  during  optimization.  The 
program  changes  only  the  leg  thick¬ 
nesses  during  the  first  descent  cycle 
when  going  from  the  input  design  to 
the  initial  design  on  the  constraint 
curve.  Therefore,  at  least  one*,  leg 
thickness  must  be  allowed  to  vary 
during  optimization  or  the  program 
will  not  be  able  to  find  an  initial 
feasible  design  point  since  it  will 
not  have  anything  to  vary. 


CARD  TYPE  VARIABLE 

2  NUMBER  CONTROLLED  ISA  ISB  ISC 


1  Side  Length  12 

2  Height  12 


3  Vertical  Leg 

Thickness 


4  5 


4  Skew  A  Leg  3  5 

Thickness 
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CARD  TYPE 

2  NUMBER 

VARIABLE 

CONTROLLED 

ISA 

5 

Skew  B  Leg 
Thickness 

3 

6 

Upper  Face 
Thickness 

Sheet 

12 

7 

Lower  Face 
Thickness 

Sheet 

12 

8 

THETA 

12 

9 

YIOFF 

12 

10 

Y20FF 

12 

11 

THT 

12 
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ISC 


L 


CARD 

TYPE  VARIABLE  FIELD 

3  Title  A(l-48) 


COMMENT 
Title  Card 
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CARD 


TYPE 

VARIABLE 

FIELD 

4 

EX 

F(l-12) 

EY 

F(13-24) 

GXY 

F(25-36) 

UXY 

F(37-48) 

COMMENT 

Elastic  modulus  in  longitudinal 
direction  of  leg,  psi. 

Elastic  modulus  in  transverse 
direction  of  leg,  psi. 

Shear  modulus  of  leg, psi. 

Poisson's  Ratio 


One  card  type  4  must  be  input 
for  each  leg  (Vertical  Skew  A, 
Skew  B)  and  for  the  upper  and 
lower  faces  if  "IFACE"  =  1. 
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CARD 

TYPE 


VARIABLE 


FIELD 


COMMENT 


5  LX  1(1-4)  Number  of  legs  in  X  face  of 

Tetra-core  "LY"  times  the  length 
of  a  tetrahedron  times  0.866  gives 
the  panel  width. 

LY  1(5-8)  Number  of  legs  in  Y  face  of 

Tetra-core  "LY"  times  the  length 
of  a  tetrahedron  side  gives  the 
panel  length. 

SIDE  F (13-18)  Length  of  a  tetrahedron  side,  in. 

HT  F (19-24)  Height  of  Tetra-Core  Panel,  in. 

THK(l)  F (25-30)  Thickness  of  Vertical  Legs,  in. 

THK(2)  F (31-36)  Thickness  of  Skew  A  Legs,  in. 

THK(3)  F (37-42)  Thickness  of  Skew  B  Legs,  in. 

THK(4)  F(43-48)  Thickness  of  Lower  Face  Sheet,  in. 

THK(5)  F (49-54)  Thickness  of  Upper  Face  Sheet,  in. 

"THK(4)"  and  "THK(5)"  are  input 
only  if  "IFACE"  =  1  on  card  typo 

1. 
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CARD 

TYPE 

VARIABLE 

FIELD 

COMMENT 

6 

THETA 

F(l-12) 

The  variables  input  on  this 
card  control  the  variation  of 
of  the  Tetra-core  element  from 

YIOFF 

F(13-24) 

a  true  tetrahedron.  See 

Figure  4  for  a  description. 

Y20FF 

F(25-36) 

THT 

F{37-48) 
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CARD 

TYPE 


VARIABLE 


FIELD 


COMMENT 


7 

TC 

F(l-12) 

Thickness/Chord  Ratio 

XC 

F(13-24) 

X/Chord  Ratio 

Card  Type  6  is  input  only  if 
"ITYP"  =  3  on  card  I,  that  is, 
only  if  an  airfoil  shape  is 
to  be  generated.  Eleven  cards 
of  type  6  must  be  input,  each 
with  a  thickness/chord  for  the 
specified  X/Chord.  The  X/Chord 
at  the  nose  of  the  airfoil  is 
input  first,  with  X/Chord  for 
the  tail  of  the  airfoil  last. 
The  intermediate  X/Chord  cards 
must  be  in  sequence  from  nose 
to  tail. 


L 


CARD 

TYPE 

VARIABLE 

FIELD 

COMMENT 

8 

STTl 

F(l-12) 

Longitudinal  tension  strain 
increment ,  in . /in . 

STT2 

F(13-24) 

Transverse  tension  strain 
increment,  in. /in. 

STCl 

F(25-36) 

Longitudinal  compression  strain 
increment ,  in . /in . 

STC2 

F(37-48) 

Transverse  compression  strain 
increment,  in. /in. 

STSS 

F(49-60) 

Shear  strain  increment,  in. /in. 

"STTl",  "STT2",  "STCl",  "STC2", 
"STSS"  are  the  strain  increments 
used  in  inputting  tension  and 
compression  stress-strain  curvea 
on  cards  of  Type  9.  The  stress- 
strain  curves  are  broken  dawn 
into  increments  of  "STTy  "STT2," 
etc.,  for  inputting  on  card  9. 
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CARD 

TYPE  VARIABLE  FIELD  COMMENT 


9 


SIGT(I)  F(l-12)  Stress-strain  curves  of  the 

I  *»  1,7  F  (13-25) ..  .material  used  in  each  leg  and 

F (61-72)  each  face  sheet  are  input  if 
F (73-80)  a  nonlinear  analysis  is 

performed  ("NLIN"  =  1,  card  1), 
Longitudinal  tension  stresses 
are  input  for  strains  of  0,  STTl, 
2STT1,  3STT1, . . .  6STT1 ,  for  example. 


The  first  card  following  card  8 
must  contain  the  longitudinal 
tension  stresses;  the  second  card 
contains  the  transverse  tension 
stresses,  etc.  The  complete 
sequence  is  shown  below. 
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CARD 

TYPE  VARIABLE  CONTENTS 


8  Strain  increments  -  "STTl,"  "STT2," 

"STCl,"  "STC2/’  "STSS.",  in. /in. 

9  SIGTl  Longitudinal  tension  stresses  in  increments 

of  strain  "STTl.",  psi. 


9  SIGT2 

9  SIGCl 

9  SIGC2 

9  PRXT 


Transverse  tension  stresses  in  incr  ^nts 
of  strain  "STT2.",  psi. 

Longitudinal  compression  stresses  in 
increments  of  strain  "STCl."/  psi. 

Transverse  compression  stresses  in 
increments  of  strain  "STC2.",  psi. 

Poisson's  ratio  wXY  in  tension  in 
increments  of  "STTl." 


9  PRXC  Poisson's  ratio  yXY  in  compression 

increments  of  "STCl." 


One  set  of  Card  8  -f  six  card  9's  must  be 
input  for  each  leg  of  the  Tetra-core.  If 
"IFACE"  *  0  on  card  1,  then  three  sets 
must  be  input-  Vertical  legs,  Skew  A  legs, 
and  Skew  B  legs.  If  "IFACE"  «  1,  then  five 
sets  must  be  inputs  three  sets  for  the  legs, 
plus  one  set  for  the  lower  face  sheet, and 
one  set  for  the  upper  face  sheet. 
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CARD 

TYPE 

VARIABLE 

FIELD 

COMMENT 

10 

FXT 

F(l-12) 

Longitudinal  tension  allowable 
stress,  psi. 

FYT 

F(13-24) 

Transverse  tension  allowable 
stress,  psi. 

FXC 

F(25-36) 

Longitudinal  compression  allow¬ 
able  stresses,  psi. 

FYC 

F(37-48) 

Transverse  compression  c.llowable 
stress,  psi. 

SSS 

F(49-60) 

Shear  Allowable  Stress,  psi. 

One  card  type  9  must  be  input 
for  each  leg.  If  "I FACE"  =  1, 
five  card  9*s  are  input,  one 
each  for  the  Vertical  legs, 

Skew  A  legs,  Skew  B  legs,  lower 
face  sheet,  upper  face  sheet. 

If  "IFACE"  =  0  only  the  first 
three  are  input. 


CARD 


TYPE 

VARIABLE 

FIELD 

COMMENT 

11 

RHO (1) 

F(l-12) 

Card  10  is  input  only  if  "lOPT" 

=  1  on  card  1. 

Density  of  material  used  in 
vertical  legs,  Ib/in . ^ . 

RHO  (2) 

F(13-24) 

Density  of  material  used  in 

Skew  A  legs,  Ib/in,^. 

RHO (3) 

F(25-36) 

Density  of  material  used  in 

Skew  B  legs,  Ib/in.-^ 

RHO  (4) 

F(37-48) 

Density  of  material  used  in 
lower  face  sheet,  Ib/in •  • 

RHO (5) 

F{49-60) 

Density  of  material  used  in 
upper  face  sheet,  Ib/in*  • 

The  density  is  used  to  calculate 
the  weight  of  the  Tetra-core  for 
use  in  optimization. 
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CARD 

TYPE 

VARIABLE 

FIELD 

COMMENT 

12 

XLOAD 

F(l-12) 

Load  on  XFACE  of  flat  plate 
Tetra-core ,  Ib/in.;  or  total  end 
load  on  cylinder  or  airfoil,  lb 

YLOAD 

F(13-24) 

Load  on  YFACE  of  flat  plate 
Tetra-core,  Ib/in. 

XYLOD 

F(25-36) 

Shear  load  on  flat  plate  Tetra- 
core,  Ib/in. 

XMOM 

F (37-48) 

Moment  on  XFACE  of  flat  plate 
Tetra-core  in. -Ib/in.;  or  total 
moment  on  cylinder  or  airfoil, 
in.- lb. 

YMOM 

F(49-60) 

Moment  on  YFACE  of  flat  plate 
Tatra-core  in. -Ib/in. 

XYMOM 

F (61-72) 

Twisting  moment  on  flat  or 
plate  Tetra-core  in. -Ib/in. 

TORQ 

F(73-80) 

Torque  on  cylinder  or  airfoil 
of  Tetra-core,  in. -Ib/in. 
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CARD 

TYPE  VARIABLE  FIELD  COMMENT 

13  XQSHR  F(l-12)  Shear  load  on  XFACE  of  flat 

plate  xetra-core, Ib/in. 

YQSHR  F (13-24)  Shear  load  on  YFACE  of  flat 

plate  Tetra-core,lb/in. 

There  is  an  upper  limit  of  10 
on  the  number  of  load  cases  that 
may  be  run. 

If  a  nonlinear  analysis  is  to  be 
made  (NLIN"  =1)  or  an  optimiza¬ 
tion  is  to  be  run  ("lOPT"  =  1) , 
only  one  load  case  is  allowed. 
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CARD 

TYPE 

14 


VARIABLE  FIELD  COMJIENT 

J  I  (1  -  4)  Card  type  14  is  input  only- if 

"ILOD"  =  1  or  2  on  card  1. 

Node  which  loads  input  on  this 
card  will  be  applied  to. 


ux 

(J) 

F 

(12-24) 

Load 

in 

X 

direction 

on 

node 

J 

UY 

(J) 

F 

(25-36) 

Load 

in 

Y 

direction 

on 

node 

J 

UZ 

(J) 

F 

(37-48) 

Load 

in 

Z 

direction 

on 

node 

J 

A  card  Type  14  is  input  for  each  node 
being  loaded.  Input  of  Card  Type  14  is 
terminated  by  a  card  with  J  =  300. 
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CARD 

TYPE 

VARIABLE 

FIELD 

COMMENT 

15 

Card  Type  15  is  input  only  if 
"ILOD"  =  2  or  3  on  Card  1. 

J 

I 

(1  -  4) 

Node  for  which  fixity  is  read 

Code  (J) 

F 

(13-24) 

Fixity  of  node  J.  Convention 
shown  below  is  used. 

Code  (J) 

Freedoms 

Fixed 

0.0 

None 

1.0 

X 

2.0 

Y 

3.0 

Z 

4.0 

X, 

Y 

5.0 

X, 

Z 

6.0 

y, 

Z 

7.0 

X, 

Y, 

Z 

88 


L 


1 


Card  Formats 

Input  to  the  Tetra-core  analysis  program  must  be  in  the  form 
of  data  cards  punched  according  to  fixed  formats. 

Integer  Fields;  I  (1-4) »  I  (26-30)/  etc. 

An  integer  must  be  right  adjusted  in  the  field. 

Unneeded  field  space  may  be  left  blank.  Blanks 
are  interpreted  as  zeros  in  the  corresponding 
column.  Decimal  points  are  not  used. 


Decimal  Fields;  F(l-12),  F(37-48),  etc. 

A  decimal  number,  punched  with  a  decimal  point, 
may  be  located  anywhere  within  the  field.  Blank 
fields  are  interpreted  as  zeros. 


ahanumeric  Fields;  A(l~27) ,  A (21-30), 


Alphanumeric  characters  may  be  anywhere  in  the 
field.  Legal  characters  for  this  field  include 
all  FORTRAN  characters. 
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LIST  OF  SUBROUTINES 


OVERLAY  TETRA  1  -  ROOT  OVERLAY 

MAIN  -  Main  program.  Calls  subroutines  GEN,  STRSS, 

SOLPAC,  DEFL. 

OUTIN  -  Used  to  read  and  write  stiffness  matrix  and 

load  vectors  on  file  4.  Called  by  STIFF,  SOLPAC, 
DEFL. 

OVERLAY  TETRA  2 

GEN  -  Subroutine  to  read  and  print  input  data.  Calls 

subroutines  NODGEN,  PLATEN,  LODGEN,  NODGAN,  PLATGA, 
LODGAN,  OPTM,  NONLIN.  Called  by  MAIN. 

CURV  -  Curve  fitting  subroutine,  used  in  generation  of 
airfoil  shape.  Called  by  NODGEN,  NODGAN. 

OPTM  -  Optimization  subroutine.  Determines  direction  of 
travel  to  obtain  minimum  weight  design.  Calls 
WEIGHT,  DAR.  Called  by  GEN. 

WEIGHT  -  Subroutine  to  calculate  weight  of  Tetra-core 
model  being  optimized.  Calls  NODGEN,  NODGAN. 
Called  by  OPTM. 

DAR  -  Subroutine  to  calculate  direction  cosines  of 

gradient  given  partial  derivatives.  Called  by  OPTM. 

NONLIN  -  Subroutine  to  increment  loads  before  each  step 
of  nonlinear  analysis.  Called  by  GEN. 

OVERLAY  TETRA  3 

NODGEN  -  Subroutine  to  generate  node  points  for  true 
Tetra-core  model.  Calls  CURV. 

OVERLAY  TETRA  4 

PLATGN  -  Subroutine  to  generate  plates  for  true  Tetra- 
core  model.  Called  by  GEN. 

OVERLAY  TETRA  5 

LODGEN  -  Subroutine  to  generate  nodal  loads  and  nodal 

fixities  for  true  Tetra-core  model.  Called  oy 
GEN. 


OVERLAY  TETRA  6 


NODGAN  -  Subroutine  to  generate  node  points  for  truncated 
Tetra-core  model.  Called  by  GEN,  WEIGHT. 

OVERLAY  TETRA  7 

PLATGA  -  Subroutine  to  generate  plates  for  truncated 
Tetra-core  model.  Called  by  GEN. 

OVERLAY  TETRA  8 

LODGAN  -  Subroutine  to  generate  nodal  loads  and  nodal 

fixities  for  truncated  Tetra-core  model.  Called 
by  GEN. 

OVERLAY  TETRA  9 

STRSS  -  Overlay  to  perforin  finite-element  stress  analysis 
on  model  generated  by  previous  overlays.  Called 
by  MAIN.  Calls  STRESS,  STIFF. 

OVERLAY  TETRA  10 

STRESS  -  Subroutine  to  renumber  Tetra-core  model  to 

account  for  midside  nodes.  Called  by  STRSS. 

NID  -  Function  to  generate  new  number  for  midpoint  node. 

OVERLAY  TETRA  11 

STIFF  -  Subroutine  to  merge  element  stiffness  matrices  into 
overall  stiffness  matrix.  Called  by  STRSS.  Calls 
TRPRD,  TRIM6,  OUTIN. 

TRPRD  -  Subroutine  to  perform  matrix  multiplication 
[A]T[b][A].  Called  by  TRIM6 . 

TRIM6  -  Subroutine  to  generate  linear  strain  triangle 

and  quadrilateral  and  constant  strain  triangular 
stiffness  and  stress  matrices.  Called  by  STIFF. 
Calls  TRPRD. 

OVERLAY  TETRA  12 

SOLPAC  -  Subroutine  for  solution  of  stiffness  matrix. 

Solves  Eqn  P  =  KA  using  Choleski  triangularization 
and  back  substitution.  Called  by  MAIN. 
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OVERLAY  TETRA  13 


DEFL  -  Subroutine  for  printing  of  deflections, 

calculation  of  plate  element  stresses.  Called  by 
MAIN.  Calls  MARGIN,  MODUL,  BUKL,  PLST. 

MARGIN  -  Subroutine  to  calculate  margin  of  safety  in  each 
plate  element  using  von  Mises-Tsai  failure 
criterion.  Called  by  DEFL. 

MODEL  -  Subroutine  to  calculate  plate  secant  moduli 

based  on  strains  within  the  plate.  Called  by  DEFL. 

BUKL  -  Subroutine  to  calculate  buckling  stress  for  typical 
plate  from  each  leg.  Called  by  DEFL. 

PLST  -  Function  to  calculate  stress  given  the  strain  in  a 
plate . 

DESCRIPTION  OF  VARIABLES  USED  IN  COMMON 
COMMON /COMN/ 

NUMNP  -  Number  of  nodal  points  in  model 
NUMEL  -  Number  of  plate  elements  in  model 
NLOD  -  Number  of  load  cases 
NLIN  -  Used  to  control  noninear  analysis 
MBAND  -  1/2  bandwidth  of  stiffness  matrix,  in  freedoms 
NBAND  -  1/2  bandwidth  of  stiffness  matrix,  in  nodes 
COMMON/MATL/ 

E  -  Moduli  of  elasticity  for  each  leg 

FXT,  FYT,  FXC,  FYC,  SSS  -  Allowable  stresses  of  material 

in  each  leg 

CMS  -  I-largin  of  safety  in  each  leg 
COMMON/LODA/ 

XLOAD,  YLOAD,  XYLOD,  XMOM,  YMOM,  XYMOM,  TORQ,  XQHR, 

YQSHR  -  input  applied  loads 
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COMMON/ONE/ 


SIDE,  HT  -  Tetrahedron  side  length  and  actual  height 
(see  Figure  4) 

THK  -  Thickness  of  each  leg 

THETA,  YIOFF,  Y20FF  -  Parameters  to  control  variation 

from  theoretical  equalateral  tetra¬ 
hedron  (see  Figure  4) 

THT  -  Theoretical  height  at  which  tetrahedron  sides  would 
meet  (see  Figure  4) 

DUM  -  Not  used 

TC,  XC  -  Thickness/chord  and  X/chord  ratios  used  to 
generate  airfoil  model 

HED  -  Title 

NX,  NY  -  Number  of  legs  on  X  and  Y  face.  Called  LX  and  LY 
on  input  sheets  (see  Figure  5) 

lOPT  -  Used  to  control  optimization 

ITYP  -  Type  of  model  to  be  run — flat  plate,  cylinder, 
airfoil 

IFACE  -  Control  of  face  sheet  generation 

IHOLE  -  Used  to  control  if  hole  will  be  placed  in  model 
or  not 

ITET  -  used  to  control  type  of  tetrahedron  generated — 
true  or  truncated 

ICNT  -  Counter  for  number  of  steps  that  have  been  run 
in  nonlinear  or  optimization  analysis 

NNEW  -  Used  to  determine  which  leg  a  given  plate  is  in 

COMMON/OPTA/ 

NSTP  -  Limiting  number  of  steps  for  nonlinear  analysis, 
limiting  number  of  cycles  for  optimization 

IPRT  -  Controls  printing  of  intermediate  printout  in 
optimization 
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IDS  -  Controls  steepest  descent  mode  in  optimization 

IDRV  -  Controls  calculation  of  derivatives  in 
optimization 

ITOT  -  Counter  for  number  of  cycles  in  optimization 
ICYC  -  Counter  for  number  of  sidesteps  in  optimization 
L^ST  -  Not  used 

IK  -  used  to  identify  variable  being  incremented  when 
calculating  derivatives  in  optimization 

WTL  -  Base  weight  before  start  of  calculation  of 
derivatives  in  optimization 

XMSL  -  Margin  of  safety  at  start  of  last  sidestep 
in  optimization 

XMS  -  Present  margin  of  safety 

WTMN  -  Minimum  weight  found  in  optimization 

TMLT  -  Step  length  used  in  sidestepping  -  optimization 

WT  -  Present  weight 

HMAX  -  Value  of  variables  associated  with  WTMN  -  optimization 

HL  -  Value  of  variables  at  start  of  calculation  of 
derivatives  -  optimization 

CMSA  -  Value  of  margins  of  safety  at  start  of  calculation 
of  derivatives  -  optimization 

SLW  -  Partial  derivatives  of  effect  of  variables  on  weight  - 
optimization 

SLOPE  -  Partial  derivatives  of  effect  of  variables  on 
margin  of  safety  -  optimization 

RHO  -  Density  of  material  in  each  leg 

PWT  -  Gradient  of  constant  weight  surface  -  optimization 

PMT  -  Composite  constraint  gradient  -  in  optimization 

ISA,  ISB,  ISC  -  Used  to  control  variables  during 
optimization 


94 


COMMON/MODU/ 

SIGTl,  SIGT2,  SIGCl,  SIGC2,  SIGSS  -  Stresses  defining 

stress-strain  curves  used  in  nonlinear  analysis 

PRXT,  PRXC  -  Poisson's  ratios  used  in  nonlinear  analysis 

STTl,  STT2 ,  STC2,  STSS  -  Strain  increments  used  to  define 
stress-strain  curves  for  nonlinear  analysis 

COMMON/LODB/ 

APPLD  -  XLOAD,  YLOAD,  etc.  for  each  load  case  are  stored 
in  APPLD 

COMMON/RAX/ 

NR4  -  Pointer  for  next  record  to  be  read  on  File  4 

NR4L0D  -  Gives  location  of  first  record  in  loads  vector 
on  File  4 

NR4STF  -  Gives  location  of  first  record  in  stiffness 
matrix  on  File  4 

NR4DIA  -  Gives  location  of  first  record  in  diagonal  on 
File  4 

NREC  -  Number  of  records  required  to  store  throe  rows  of 
stiffness  matrix  on  File  4 

NRECST  -  Number  of  records  in  stiffness  matrix 

NRECLD  -  Number  of  records  in  loads  vector 

NRl  -  Pointer  for  next  record  to  be  read  on  File  1 

NR2  -  Pointer  for  next  record  to  be  read  on  File  2 

DESCRIPTION  OF  VARIABLES  USED  IN  SUBROUTINES 

SUBROUTINES  -  GEN,  NODGEN,  PLATGN,  LODGEN ,  NODGAN ,  PLATGA, 
LODGAN,  WEIGHT 

X,  y,  Z  -  Nodal  coordinates 

SLOPA  -  Slope  of  surface  at  each  node--used  for  airfoil 
NA,  NB,  NC,  ND  -  Nodes  at  vertices  of  plate  element 
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CODE  -  Fixity  for  each  node 
UX,  UY,  UZ  -  Loads  on  each  node 
SUBROUTINE  -  OPTM 

SLP  -  Intermediate  step  in  calculation  of  composite 
constraint  gradient 

PMS  -  Gradient  to  each  margin  of  safety  surface 

il  -  Vector  of  variables  <o  be  optimized 

SUBROUTINE  -  STRESS,  STIFF,  TRIM6 

X,Y,Z,UX,UY,UZ,CODE  -  See  subroutines  NODGEN ,  PLATGEN ,  etc. 

IDO  -  Node  numbers  as  renumbered  to  include  mid-side  nodes 

ID  -  Node  numbers  in  oiiginal  ordoj  ,  not  including 
midside  nodes 

NI  -  Node  number  of  nev;  mid-side  node 

IN  -  Senuenco  of  new  mid-side  nt^de 

ML  -  Added  mid-side  nodes  in  plate  element 

LM  -  Original  vertex  nodes  inpiate  eJoment  -  in  original 
nodal  n'ornbering  sequence 

Mil  -  Renumbered  mid-side  and  vertex  nodes  in  plate  element 

DUMl , DLM2 , DUM3 , DUM4  -  Dummy  variables  used  in  renumbering 

of  load  vector  and  nodal  fixities  to 
include  the  effect  of  added  midside 
nodes 

S  -  Element  stiffness  matrix 

ETR,ELGT,EPRM,  G  -  Plate  elasticity  matrix  -  ETR  =  EY/(1- 

II  ),  ELGT  =  Ex/(l-y  n  ),  EPRM  = 
xy  yx  xy  yx 

Ey,;xy/ ( 1-iJxyMyx)  ,  G  =  GXY 


T  -  Plate  thickness 

ST  -  Plate  stress  matrix 

IMTL  -  Indicates  leg  that  plate  is  in 
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A  -  Three  rows  of  merged  stiffness  matrix 

DIAG  -  Diagonal  of  merged  stiffness  matrix 

PHYX,  PHYY  -  Used  in  computing  element  stiffness  matrix 

ALAM  -  Element  transformation  matrix.  Transforms 

stiffnesses  from  local  to  global  coordinates 

AL,  AL2  -  Direction  cosines  of  element 

QMAT  -  Elasticity  matrix 

AST,  AK21,  DUM,  STA,  STR,  STRA  -  Dummy  matrices  used  in 

calculating  stiffness  and 
stress  matrices 

ZETA  -  Area  coordinates  of  location  in  triangle  for 
which  stresses  are  calculated 

IK  -  Order  in  which  triangle  stiffnesses  are  added  into 
quadrilateral  stiffness  matrix 

SUBROUTINE  SOLPAC 

DIAG  -  Diagonal  of  stiffness  matrix 

AI  -  Three  rows  of  stiffness  matrix 

AJ  -  Three  columns  of  stiffness  matrix 

B  -  Initial  loads  vector,  deflections  after  solution 
is  complete 

X  -  Loads  vector  after  reduction 
SUBROUTINES  DEFL,  BUKL,  MODUL 

P  -  Displacements  for  nodes  connecting  to  an  element 
LM,  MM,  ST,  IDO  -  See  subroutine  STRESS 
STRESS  -  Element  stresses 

FBKL  -  Buckling  stress  of  an  element  in  each  leg 
STRAIN  -  Element  strains 
B  -  Nodal  deflections 
D  -  Plate  stiffnesses 
XL  -  Width  of  plate 
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Merged  stiffness  matrix 

Triangularized  stiffness  matrix 

Matrix  relating  plate  strains  to  plate  nodal 
deflections 

Plate  elasticity  matrix 

Plate  bending  stiffnesses 

Plate  elastic  modulus  in  fiber  direction 

Plate  elastic  modulus  transverse  to  fiber 
direction 

Plate  tension  allowable  stress  in  fiber 
direction 

Plate  compression  allowable  stress  in  fiber 
direction 

Plate  tension  allowable  stress  transverse 
to  fiber  direction 

Plate  compression  allowable  stress  transverse 
to  fiber  direction 

Plate  shear  modulus 

Number  of  Vertical  legs  on  X  face 

Number  of  Skew  B  legs  on  Y  face 

Plate  buckling  load  in  X  direction  (Ib/in.) 

Vertex  nodal  load  for  uniform  load  in  X 
direction  ’ 

Vertex  nodal  load  for  uniform  load  in  Y 
direction 

Vertex  nodal  load  for  uniform  load  in  XY 
direction 

Vertex  nodal  load  for  uniform  moment  in  X 
direction 
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MY 

-  Vertex  nodal 

direction 

^MXY 

-  Vertex  nodal 

direction 

PZ 

-  Vertex  nodal 

lead  for  uniform  moment  in  Y 

load  for  uniform  moment  in  XY 

load  in  Z  direction 
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APPENDIX  I 

ADDITIONAL  GENERATION  EXAMPLES 


Additional  examples  of  flat  plate  model  generation  are  included 
in  this  appendix  to  clarify  the  examples  given  in  the  body  of 
the  report.  The  nodal  numbering  systems  used  for  a  6  by  6  and 
a  10  X  6  true  Tetra-core  plate  are  shown  in  Figures  29  and  30, 
respectively.  The  nodal  numbering  systems  used  for  a  4  x  6  and 
an  8  X  4  truncated  Tetra-core  plate  are  shown  in  Figures  31  and 
32.  Plate  numbering  systems  used  for  the  Vertical,  Skew  A, and 
Skew  B  legs  in  the  truncated  Tetra-core  flat  plate  are  shown 
in  Figures  33,  34,  and  35.  The  addition  plates  added  to  Skew  A 
and  Skew  B  legs  to  complete  the  cylinder  model  are  shown  in 
Figures  36  and  37.  The  nodal  loads  generated  for  a  true  Tetra- 
core  flat  plate  model  for  each  load  type  are  shcvm  in  Figures 
38,  39,  40,  41,  42,  43,  44,  and  45.  The  nodal  loads  generated 
for  a  truncated  Tetra-core  flat  plate  model  are  shown  in  Figures 
46,  47,  48,  49,and  50. 


V 
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•  Nodes  Connected  by  Dark  Lines  Are  on  Upper  Surface  of  Plate 

•  Nodes  Connected  by  Light  Lines  Are  on  Lower  Surface  of  Plate 


Figure  29.  Nodal  Numbering  System,  True  Tetra-Core  Flat  Plata, 
*  6,  Ly  ^  6. 
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Figure  30.  Nodil  Numbering  Syttem,  True  Teire-Core  net  Plete.  Lj(  •  10.  L 


Z  •  3 
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Figure  31.  Nodal  Numbering  System,  Truncated  Tetra-Core  Flat  Plate, 
Lx  =  4,  Ly  =  6. 
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Figura  32.  Nodal  Numbaring  Systam,  Truncatad  Tatra-Cora  Flat  Plata, 
Lv  *  8,  Lw  ^ 


Figure  33b  Verticil  Leg  Plele  Numbering  SyitHn  Trunceled  Tetre-Core. 

Lx  -  4,  Ly  -  4. 
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First  Skew  A  Plate 


X 


Figur*  34.  Skew  A  Leg  Plate  Numbering  System  Truncated  Tatra-Core, 
Lx  =  4,  Ly  =  4. 
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Fifurt  3S.  Skew  B  LegPleee  Numbering  Syetwn  Truncated  Tctra-Core, 

Lx  *  4,  Ly  -  4. 
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Fir«  Plata  in 
SkMv  A  Lags 


Skew  A  to  "Zip  Up" 
Side  of  Cylinder 
or  Airfoil 


F»gur«  36. 


A  Pim  NumlMring  Syitwn,  True  TffrMJofg 

‘■X  *  6,  Ly  »  6. 


FlitPIm. 


189 


Last  Plate  in  Skew  B  Legs 


Skew  B  to  "Zip  Up' 
Side  of  Cylinder  or 
Airfoil 


Sk«w  B  Ptali  Numbering  Sy«lMn,  True  Tetre-Core  Flat  Plate, 
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Fi^ir*  3&  Load  Appliad  to  Tnjt  Tatra-Cort  Flat  Plata,  L 


Fi^a  39.  Load  Ny  Applied  to  True  T«fr»-Cor«  Flat  Plata, 


Figure  40.  Ir>-Piarte  Shear  Load  Appliad  to  Trua  Tatra-Cort  Flat  Plata. 
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Figure  ^  Moment  My  Applied  to  True  Tetra-Core  Flat  Plate, 


Fitur«43b  Twitting  Moment  M^y  Applied  to  Tru«  T«tr»-Cor«  Flat  Plata, 


PMX  -  MooMnt  Lotd*  to  Baionco  PZ  Loads 


Varticil  ShMr  Load  XQSHR  Appliad  to  Tatra-Cora  Flat  Plata,  Lx 


fZ  -  Vtrtical  ShMr  Nodal  Load*  PMY  -  Momant  Loads 


•  =  ♦Z  Load 
□  =  -Z  Load 


Fipurt45b  Varticil  ShMr  Lodd  YQSHR  AppUsd  to  T«Cra-Cor«  Flat  Platt,  -  4,  Ly 
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Figure  4&  Load  Applied  to  Truncated  Tetra-Core  Flat  Plate,  ■  4,  Ly 
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Figurt47.  Loid  Ny  Appliad  to  Truncatad  TttnhCor*  Flat  Plata.  Ly  *  4.  Ly  4. 
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Fipm  41.  In-Plm  Star  Load  N^y  Appliad  to  Truncatad  Tair»<€ora  Flat  PMa, 
L)(  ■  4,  Ly  ■  4. 
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Fi|iir«4§.  MofiMntM^  AppliidtoTninGatadT«tr»-CcM  FlatPlaii,  L)(  ■  4,  Ly  ■  4. 
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Figure  50.  Moment  My  Applied  to  Trunceted  Tetre-Core  Flet  Plete,  L 


APPENDIX  II 

EXAMPLE  OF  COMPUTER  PROGRAM  OUTPUTS 


An  excunple  of  the  output  from  an  optimization  run  is  given 
(Figure  51)  .  Stresses  and  deflections  for  the  input  configura¬ 
tion  are  printed.  Then  the  optimization  steps  are  printed. 
Finally,  stresses  and  deflections  for  the  optimized  section  are 
printed.  Element  stresses  are  printed  out  of  sequence  in  this 
example,  although  they  are  in  sequence  in  the  latest  version 
of  the  progreun. 
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loads  opt  Plastic  ttpc  pace  hole  steps  tetra 
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Fifuro  51.  Output  frofn  Optimization  Run. 
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